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Abstract. 

In the first part of this review, experimental results which lead to the discovery 
of two-proton radioactivity are examined. Beyond two-proton emission from nuclear 
ground states, we also discuss experimental studies of two-proton emission from excited 
states populated either by nuclear (3 decay or by inelastic reactions. In the second 
part, we review the modern theory of two-proton radioactivity. An outlook to future 
experimental studies and theoretical developments will conclude this review. 
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1. Introduction 

Atomic nuclei are made of two distinct particles, the protons and the neutrons. These 
nucleons constitute more than 99.95% of the mass of an atom. In order to form a stable 
atomic nucleus, a subtle equilibrium between the number of protons and neutrons has 
to be respected. This condition is fulfilled for 259 different combinations of protons and 
neutrons. These nuclei can be found on Earth. In addition, 26 nuclei form a quasi stable 
configuration, i.e. they decay with a half-life comparable or longer than the age of the 
Earth and are therefore still present on Earth. 

These stable or quasi-stable nuclei are characterised by an approximately equal 
number of neutrons and protons for light nuclei, whereas going to heavier nuclei the 
neutron excess increases to reach N-Z=44 for the heaviest stable nucleus ^°^Pb. In 
addition to these 285 stable or quasi-stable nuclei, some 4000-6000 unstable nuclei are 
predicted to exist, the exact number depending on the theoretical model used. Close 
to 2500 nuclei have been observed already. These unstable radioactive nuclei can be 
classified into seven categories depending on their decay modes: i) a emitters (about 
375) which decay by ejecting a ^He nucleus, ii) emitters (1040) which transform a 
proton into a neutron plus a positron and a neutrino, or by capturing an electron from 
their atomic shells and emitting a neutrino, iii) (3" emitters (1020) which transform 
a neutron into a proton by emitting an electron and an anti-neutrino, iv) fissioning 
nuclei (30) which split up into two roughly equal pieces, v) one-proton (Ip) emitters 
(about 25), vi) two- proton (2p) emitters (2) and, finally, (vii) exotic-cluster emitters 
(about 25 different nuclei emitting fragments of C, O, F, Ne, Mg, and Si). These exotic 
radioactivities with masses from 14 to 32 may be considered as intermediate between 
a decay and fission. However, it should be mentioned that this decay mode is never 
dominant. 

Although some radioactive nuclei can decay by two or three of these different decay 
modes, the main decay mode is dictated by the proton-to-neutron ratio. For example, 
nuclei with a modest proton excess decay by disintegration. If the proton excess 
increases, the nuclear forces can no longer bind all protons and one reaches the proton 
drip line where the spontaneous emission of protons takes place from the ground state. 
Heavier neutron-deficient isotopes will preferentially decay by a emission, and for the 
heaviest nuclei fission is the dominant decay mode. For some of the heaviest nuclei, 
exotic-cluster emission is a small decay branch. 

The study of these different decay types started with the discovery of radioactivity 
by H. Becquerel in 1896. Following the work of Pierre and Marie Curie, E. Rutherford in 
1899 classified the observed radioactivity into two different decay modes (a and (3 rays) 
according to the detected particle in the decay. P. Villard discovered 7 radiation in 1900 
as a third radioactive decay mode. Finally, the nuclear fission process was observed for 
the first time by O. Hahn and F. Strassmann and correctly explained by L. Meitner in 
1938. 

These "classical" decay modes have helped largely to understand the structure of 
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the atomic nucleus and the forces acting in the nucleus. They have also found many 
applications which range from solid-state physics to astrophysics and medicine. With 
the discovery of a large number of new radioactive isotopes the understanding of nuclear 
structure advanced as well. 

Beginning of the 1960's, Goldanskii [1], Zel'dovich [2j, and Karnaukhov [3J proposed 
new types of proton radioactivity in very proton-rich nuclei. For nuclei with an odd 
number of protons, one-proton radioactivity was predicted. This decay mode was indeed 
observed at the beginning of the 1980's in experiments at GSI, Darmstadt [HE]. Today 
about 25 ground-state one-proton emitters are known [6]. Proton emission was also 
observed in many cases from long-lived excited states, the first being ^^'"Co [ZIIHIIS]. The 
study of one-proton radioactivity allowed testing the nuclear mass surface, to determine 
the sequence of single-particle levels and the detailed structure of the wave function 
of the emitted proton ("j content" of the wave function), and to investigate nuclear 
deformation beyond the proton drip line. 

Zel'dovich [2] was probably the first to mention the possibility that nuclei may emit 
a pair of protons. However, it was Goldanskii [U |10l HIl |T2l |T3] and Janecke [H] who 
first tried to determine candidates for 2p radioactivity. Many others followed and the 
latest predictions [l5l[T6l[n] determined that ^^Ti, ^^Cr, ^^Fe, and ^^'^^Ni should be the 
best candidates to discover this new decay mode. Goldanskii [1] also coined the name 
"two-proton radioactivity" and it was Galitsky and Cheltsov [H] who proposed a first 
theoretical attempt to describe the process of 2p emission. More than 40 years after its 
theoretical proposal, 2p radioactivity was discovered [IHl [20] in the decay of "^^Fe. In 
later experiments, ^^Zn [21] and most likely ^^Ni [22] were also shown to decay by two- 
proton radioactivity. The observation of two-proton emission from a long-lived nuclear 
ground state was preceded by the observation of two-proton emission from very short- 
lived nuclear ground state, e.g. in the cases of ^Be j23j and ^^O [21], and by the emission 
of two protons from excited states populated by /? decay |25] or inelastic reactions [26] . 
A particular case is the emission of two protons from an isomer in ^'^Ag'" [27] . 

The 2p emission process leads from an initial unbound state to a final state of 
separated fragments interacting by the Coulomb force. This transition to the final state 
can be realized in many ways, such as the sequential two-body decay via an intermediate 
resonance, the virtual sequential two-body decay via the correlated continuum of an 
intermediate nucleus, the direct decay into the three-particle continuum or its particular 
diproton ('^He') limit of two sequential two-body emissions. 

Different types of theoretical approaches exist which focus on the three-body decay 
width: the extended R- matrix theory [28l |29] , the real-energy continuum shell model 
(the so-called Shell Model Embedded in the Continuum (SMEC) [301 [3ll [321 [33l [M] ) 
describing sequential two-body emission, virtual sequential two-body emission and 
diproton emission via two sequential two-body emissions [33l |3l], the Gamow Shell 
Model (GSM) with no separation between two- and three-body decays from the multi- 
particle continuum [35], the three-body models with outgoing flux describing direct 
decay into the continuum [36] , the three-body models combined with complex coordinate 
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scaling [23, or the Faddeev equations combined with either outgoing flux or complex 
scaling [38] . 

In the 2p decay, the three-body decay width provides global information with 
weighted contributions of different decay paths leading from the initial to the final 
state. Hence, the value of the width by itself does not reveal the dominant transition 
mechanism. In the theoretical analysis however, one can switch off certain decay paths 
and look for the three-body partial decay width resulting from a particular decay path. 
This strategy is used in the SMEC analysis of the experimental results. Obviously, the 
conclusions of such partial studies of decay mechanisms are not unambiguous. 

In the present review, we will summarize the state of the art of experimental and 
theoretical studies on two-proton radioactivity. Previous reviews, written before the 
actual discovery of ground-state two-proton radioactivity, can be found in [39l |40] . 

2. Basic concepts for one- and two-proton radioactivity 

Protons are charged particles, therefore they are sensitive to the charge of other protons 
which create a Coulomb barrier. This barrier prevents protons from quickly leaving the 
atomic nucleus even if they are unbound. The tunnelling probability depends on the 
available energy and the height of the Coulomb barrier, which in turn depends on the 
nuclear charge Z (number of protons) (see Fig. [1]). The barrier penetration can give rise 
to measurable half-lives, if a certain balance between the available decay energy and the 
barrier height is respected. Fig. [2] shows in a simple model the relation between barrier 
penetration half-life and decay energy for different nuclear charges Z. In general, the 
higher is the available energy, the shorter is the tunnelling time. In turn, for higher Z 
more energy is needed for the same tunnelling time. 



Figure 1. Due to their charge the protons have to tunnel through the Coulomb 
barrier generated by other protons. For one-proton emission (left), the tunnelling 
depends mainly on the barrier height. For two-proton emission (right) the correlation 
between the two protons most likely influences the tunnelling process. Figure courtesy 
of J. Giovinazzo. 
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The delay associated with the tunnelling process allows for the observation of Ip 
and 2p radioactivity. Even if protons are unbound by e.g. 1 MeV, the tunnelling of 
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the combined Coulomb and centrifugal barriers is not instantaneous, i.e. the nuclear 
decay is delayed by a measurable amount of time. However, due to experimental 
constraints, mainly linked to the techniques used to study Ip or 2p radioactivity and 
due to the competition with f3~^ decay, observation limits exist. The lower half-life limit 
of about 1 /is comes from the fact that often the observation of the Ip or 2p emitter is 
accomplished with the same detection setup which is used to detect the decay of these 
nuclei. This means that a decay signal of typically 1 MeV has to be observed a very short 
time after an identification signal (i.e. an implantation signal in a silicon detector) of 
several hundred MeV. This observation limit has been reached in Ip-emitter studies (see 
Sect. 3). The upper limit (see Fig. [2]) strongly depends on the structure of the decaying 
nucleus which governs the /3-decay half-life. It can range from a few milliseconds to 
a few seconds. To observe charged-particle emission, the barrier-penetration half-life 
should be comparable or shorter than the /3-decay half-life. 




0.2 Oi. 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 

decay energy (MeV) 



Figure 2. Barrier penetration lialf- lives for a proton as a function of the nuclear charge 
and the decay energy. The half-lives are calculated from Coulomb wave functions using 
the Wigner single-particle width (see e.g. [E]). The horizontal line gives the lower 
detection time limit, whereas the hatched area gives typical /3-decay half-lives. 

In the Ip decay, there are only two particles in the final state and the decay is a 
simple back-to-back decay, where the energy is shared between the two partners, the 
heavy recoil and the emitted proton, according to energy and momentum conservation. 

In the case of 2p decay, the situation is more complicated (see Fig. [1]). The decay 
characteristics depend sensitively on the decay pattern itself. Two schematical pictures 
are usually used to represent possible limiting cases: i) the three-body decay and ii) the 
diproton decay. In the former case, the two protons do not have any correlation beyond 
the phase-space constraints, which means that only energy and momentum conservation 
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have to be respected. Such a decay pattern yields an isotropic angular distribution of 
protons which share the total decay energy, with individual proton energies ranging 
from zero to the total decay energy. A more realistic discussion of the tunnelling process 
changes this picture because the barrier penetration strongly favours an emission of two 
protons with similar energies. In the latter case of diproton emission, one assumes that a 
pre-formed '^He'-cluster penetrates the Coulomb barrier and decays outside the barrier. 
The decay half-life depends sensibly on the ^He resonance energy. The diproton decay 
corresponds to two sequential binary decays for which the kinematics is rather easy. 

These two limits of 2p decay are not very realistic, but they are easy to grasp and 
give at least a schematic idea about the 2p emission. A theoretical description of 2p 
radioactivity (see Sect. 8) is more sophisticated, although none of the models provide 
an unconstrained description of the nuclear structure and the decay dynamics involved 
in the 2p decay. 

Up to now, only total decay energies and half-lives have been measured in 2p 
radioactivity experiments. A deeper insight into the decay mechanism can be provided 
by the measurements of individual proton energies and the proton-proton angle in the 
center-of-mass. To access these observables, specific detectors have to be developed (see 
Sect. 6). 

3. The discovery of one-proton radioactivity 

Soon after the prediction of proton emission from atomic nuclei, the search for these new 
types of radioactivity started. The first exotic decay mode discovered was the /3-delayed 
emission of one proton in the decay of ^^Si [1^ . A similar type of proton emission from 
an excited state was observed in the search for /5-delayed proton emission from ^^Ni. 
Instead of populating ^^Ni, the authors [3, E] observed the population and decay of 
^^Co™. This isomer {Ti/2 = 247 ms) decays with a branching ratio of 1.5% by emission 
of a 1.59(3) MeV proton. 

The first case of ground-state one-proton radioactivity was reported by Hofmann 
et al. [1] and by Klepper et al. [5]. Both experiments were conducted at GSI, the first 
one at the velocity filter SHIP using a fusion-evaporation reaction with a ^^Ni beam and 
a ^^Ru target to produce the proton emitter ^^^Lu, the second one at the on-line mass 
separator with the reaction ^^Ni + ^^Mo to synthesize ^^^Tm. 

One-proton radioactivity studies have developed soon into a powerful tool to 
investigate nuclear structure close to and beyond the proton drip line |l2l HSl E] with 
now more than thirty different nuclei known to emit protons from their ground states or 
low-lying isomeric states. These studies have given valuable information on the sequence 
of single-particle levels and their energies in the vicinity of the proton drip-line, on the 
j content of the nuclear wave function, on the nuclear deformation, and they allow to 
test predictions of nuclear mass models in this region. Most of this information can be 
obtained only by means of Ip radioactivity studies. 

More recently proton radioactivity was also used to identify nuclei and to tag 
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events. In experiments e.g. at the Fragment Mass Separator of the Argonne National 
Laboratory, the observation of a decay proton is used as a trigger to observe 7 rays 
emitted in the formation of the proton emitter. Therefore, not only the decay of these 
exotic nuclei can be studied, but also their nuclear structure in terms of high-spin states 
and their decay via yrast cascades [B]. 

The relatively large number of known ground-state proton emitters (for almost all 
odd-Z elements between Z=50 and Z=83, at least one proton-radioactive nucleus is 
experimentally observed) is due to the fact that the Ip drip line, i.e. the limit where 
odd-Z nuclei can no longer bind all protons, is much closer to the valley of stability 
than the two-proton drip line. As an example, we mention that at Z=50 the T^ = — 1 
isotope ^'^Sn is most likely the first 2p unbound nucleus. For Z=51, ^°^Sb with Tz=l/2 
is probably proton unbound. This means that in this case the 2p emitter is three 
mass units more exotic than the Ip emitter, which leads to orders of magnitude lower 
production cross sections. 

4. Two-proton emission from very short-lived nuclear ground states 

The first experimental investigation of a possible ground-state 2p emitter was performed 
by Karnaukhov and Lu Hsi-T'ink [15]. These authors tried to produce ^^Ne by 
bombarding a nickel target with a 150 MeV ^''Ne beam from the JINR, Dubna 300 cm 
cyclotron. The non-observation of any 2p event led the authors to conclude that either 
the production cross section is much smaller than assumed (smaller than 10~^ b) or that 
the half-life of ^^Ne is shorter than 10~*s. As early as in 1963, Janecke predicted its half- 
life to be of the order of lO^^^s [36] and recently its half-life was given as 9xl0^^^s j47j . 

4.I. The decay o/^Be 

If we omit early measurements which used vr-induced reactions in bubble or spark 
chambers [HI HH], and mass excess measurements performed by Ajzenberg-Selove et 
al. [50] , the first successful study of a radioactive decay with an emission of two protons 
was performed in experiments with ^Be [51j. In these studies, ^Be was produced in 
a ^Li(^IIe,t)^Be reaction. The triton as well as the a particle and the protons from 
the decay of ^Be were observed. The ground state of ^Be is known to have a life-time 
of the order of 10~^^s. Therefore, this nucleus can only be identified via its decay 
products. The authors observed an enhancement of low-energy a particles compared to 
a simple phase-space distribution which could not be explained by means of final-state 
interactions. However, no conclusions on the decay mechanism could be drawn from 
these observations. 

A much more detailed study of the ^Be decay was performed over several years at 
the Kurchatov Institute, Moscow (see e.g. [23l[52]). These authors measured energy and 
angular distributions of the emitted particles and interpreted their data using Jacobi 
coordinates. The most important conclusion of this analysis is that the spectra can be 
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interpreted as a "democratic" three-body decay, i.e. the decay where none of the two- 
body subsystems of this three-body final state has a width which is small in comparison 
with the transition energy between the initial nucleus and the subsystem. 

4.2. The decay of^^O 

The next step in the investigation of short-lived 2p ground-state emitters was the study 
of ^^O at the National Superconducting Cyclotron Laboratory (NSCL) of Michigan 
State University, East Lansing, USA [23]. In this experiment, ^^O was produced by 
one- neutron stripping from a ^^O secondary beam. The complete reaction kinematics 
was reconstructed by measuring the momentum vectors of two protons in coincidence 
with the momentum of the heavy ^'^C residue. From these measurements, the authors 
could produce an excitation energy spectrum of ^^O, determine the proton-proton energy 
difference and the proton-proton angle in the ^^O center of mass. 
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Figure 3. Experimental spectra for the decay of ^^O and level scheme. The 
experimental data are compared with Monte-Carlo simulations using the two extreme 
pictures, i.e. sequential (full line) or diproton (dashed line) decay. Quantitative 
agreement is obtained within the sequential decay picture. The decay scheme evidences 
that the sequential path is energetically open [211 ■ Figures re- used with permission 
from APS. 



These experimental data were interpreted by means of Monte-Carlo simulations, 
which employed the R-matrix theory [5l| to describe the 2p decay in two extreme 
pictures: diproton emission with a ^He resonance energy of 50 keV and three-body 
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decay which is only restricted by the available phase space. Fig. [3] summarizes these 
results. A comparison of the experimental results and the simulations clearly favours 
an uncorrelated three-body emission or a sequential decay, which may have roughly 
the same characteristics. At the time of this work, these conclusions were somewhat 
surprising, as it was believed that the intermediate Ip daughter state in ^^N was 
inaccessible due to energy conservation. In addition, a very large ^^O decay width 
was needed to explain the experimental data. However, later measurements [551 IMl 157] 
revealed that the ^^N ground state lies lower in mass and opens therefore the way for a 
sequential decay. With these new results, a consistent description of the decay of ^^O 
could be provided [53] . 

4.3. The decay o/^^Ne and ^^Mg 

The decay of ^^Mg was recently studied in an experiment at the fragment separator 
FRS of GSI [58] • In this experiment, ^^Mg was produced by a one-neutron stripping 
reaction at the FRS intermediate focal plane of ^''Mg, produced itself at the FRS target 
in a fragmentation reaction of the ^^Mg primary beam. The protons emitted in the 
decay of ^^Mg were detected and tracked by a set of three DSSSD, whereas the heavy 
decay product, ^^Ne, was identified and its momentum measured in the second half of 
the FRS. 

The measured data were compared to Monte Carlo simulations and allowed to 
extract the half-life of ^^Mg of T1/2 = 4.0(15) ps and the associated two-proton Q- value: 
Q2p = 0.75(5) MeV. This value is much smaller than standard predictions like the one of 
the Atomic Mass Evaluation [17| (Q2p = 2.00(25) MeV), the Garvey-Kelson relation [59] 
{Q2p = 1.02(2) MeV), or the prescription from Antony et al. [60] (Qzp = 1.19(3) MeV). 
However, the correlation between the Q value and the decay half-life is well reproduced 
within the three-body model of Grigorenko and co-workers |36j . 

In the same experiment [58], the decay of ^^Ne by two-proton emission was also 
observed and the known decay Q value [61] of 1.41(2) MeV was used to test the 
experimental procedure. 

4.4- Conclusions and future studies with short-lived ground-state two-proton emitters 

For the systems ^Be, ^^O, and ^^Ne, the decaying states have large widths so that the 2p 
emitter state and the Ip daughter state overlap. Therefore, the sequential decay channel 
is opened and the decay will most likely take place sequentially. The situation in the 
case of ^^Mg is less clear, as the structure and the mass of the one-proton daughter ^^Na 
is not known. However, with increasing nuclear charge, the Coulomb barrier becomes 
stronger and the wave function is more and more confined in the nuclear interior which 
yields narrower states. 

Therefore, it might well be that in nuclei like ^^Si, ^^S, or ^°Ar the Coulomb barrier 
is already strong enough so that, with a possibly small decay energy, narrow levels avoid 
overlapping and the sequential decay path is, at least in a simple picture, not open. 
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5. Observation of two-proton emission from excited states 

Generally speaking, the 2p emission can take place as a sequential-in-time process or a 
simultaneous process. In its standard form, the sequential process happens when levels 
in the Ip daughter nucleus, both resonant and non-resonant ones, are accessible for Ip 
emission. This situation may be encountered when the emission takes place from highly 
excited levels of the 2p emitter. These levels can be fed either by (3 decay or by a nuclear 
reactions. 

A clear signal for a sequential emission can be obtained when the intermediate level 
through which the emission passes is identified. This is the case, for example, when well 
defined one-proton groups sum up to a fixed 2p energy. The one-proton groups give 
then the energy difference between two well-defined nuclear levels. 

5.1. (3-delayed two-proton emission 

The possibility of /3-delayed 2p (/32p) emission was first discussed by Janecke pJ^J. He 
concluded that there should be no "/^-delayed diproton emitters" with Z<14. This decay 
mode was also studied by Goldanskii [62] who proposed possible candidates. 

The first experimental observation of /9-delayed 2p emission was reported by Cable 
et al. in 1983 |25]. Their experiments, using a ^He induced reaction on a magnesium 
target and a helium-jet technique, led to the first observation of ^^Al by means of its 
/5-delayed Ip emission from the isobaric analogue state (IAS) in ^^Mg to the ground 
and first excited states of ^^Na [63]. The /32p decay mode was subsequently observed 
with two silicon detector telescopes which allowed to measure the energy of the two 
individual protons with good resolution [25] . 

From the observation of well defined one-proton groups. Cable et al. concluded that 
this decay is a sequential mode via intermediate states in the Ip daughter nucleus ^^Na. 
This conclusion was confirmed in a further experimental study [M] where it turned out 
that the angular distribution of the two protons is compatible with an isotropic emission, 
although a small angular correlation (less than 15%) could not be excluded. 

Meanwhile more /32p emitters have been identified. Shortly after the discovery of 
/32p decay from ^^Al [25], the same authors observed /32p emission from ^^P [65] and 
from ^^Ca [66]. All in all, nine /32p emitters have been observed: ^^Al [25], ^^Si [67], 

p], [68], 31 Ar [69], ^^Ca [66], ^^Ti [70], ^^Cr [THUS], and ^^Ni [73]. Most of these 
studies were characterized by rather low statistics inherent to most investigations of very 
exotic nuclei and hence no detailed search for angular correlations could be performed. 
In addition, most of the observed decays proceed via the IAS in the /3-decay daughter 
nucleus. The 2p emission from this state is forbidden by the isospin conservation law 
and can only take place via a small isospin impurity of the IAS. 

The first high-statistics study of a /32p emitter was carried out at ISOLDE in 
CERN [711 1751 1761 177] . For the first time, a rather high-efficiency, high-granularity 
setup which allows a high-resolution measurement in energy and angle was used. This 
experiment, which studied the decay of ^^Ar, demonstrated the occurrence of /32p decay 
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branches via states other than the IAS. Indeed, many Gamow- Teller fed states could 
be shown to decay by 2p emission. The main interest of these decays is that they are 
isospin allowed and cover a wider range of spins for the 2p emitter states, the possible 
intermediate states and the final states. 

Detailed studies of the proton-proton angular correlation for different 2p emission 
branches did not yield evidence for any angular correlation. In fact, the results were 
compatible with sequential emission via intermediate states in the Ip daughter nucleus. 
Although there was some activity detected with the two protons sharing the available 
decay energy equally, the statistics of these events was too low and the angular coverage 
of the setup was not uniform enough to draw any conclusion on the possible observation 
of events, where the protons were correlated in energy and angle. Therefore, despite 
a factor of 10-100 higher statistics as compared to earlier experiments, no angular and 
energy correlation, a possible signature of a strong proton-proton correlation, could be 
evidenced. 

This observation is not really astonishing. According to Brown [78], who studied 
the competition between Ip emission and correlated 2p (diproton) emission from the 
IAS in 22Mg, the weight of the direct 2p branch is expected to be at best of the order 
of a few percent and therefore rather difficult to observe experimentally. New attempts 
to search for such a correlated 2p emission in '^^Ar are on the way [79]. These studies 
take profit from a much more uniform angular coverage. 

5.2. Two-proton emission from excited states populated in nuclear reactions 

Besides populating excited 2p emitter states by (3 decay, these nuclear states can also 
be fed by nuclear reactions like pick-up, transfer, or fragmentation. In this section, we 
will describe experiments which used this kind of interaction to study 2p emission from 
excited states. 

5.2.1. Two-proton emission from ^^O In an experiment performed at Louvain-La- 
Neuve, Bain et al. [26] bombarded a (CH2)n target with a ^^N radioactive beam 
populating highly-excited states in ^^O. The 7.77 MeV state in ^^O, which is strongly 
populated in 2p transfer reactions [80j, is expected to have a significant admixture of a 
2p configuration outside of the ^^C core. In addition, this state is 2p unbound by about 
1.2 MeV and, unlike states populated by a super- allowed [3 decay, the 2p decay to the 
ground state of ^^C is isospin allowed. Therefore, this resonance was believed to have 
a small 2p decay branch, either directly to the ^^C ground state or sequentially via a 
narrow state in ^^N (see Fig. H]). 

The two protons from the decay of the 7.77 MeV resonance were detected with 
the LEDA device [26]. The data showed clearly a resonant structure when the beam 
energy was 'on- resonance'. The authors determined a 0.16(3)% 2p branch for the decay 
of this resonance. The individual energies of the two protons could be explained by 
a sequential decay pattern via the 2.37 MeV state in ^^N. Monte-Carlo simulations 
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Figure 4. Decay scheme of the 7.77 MeV resonance in ^''O by two-proton emission. 
A direct decay to the ground state of ^^C is possible, but the decay was found to be 
predominantly sequential via an intermediate state at 2.37 MeV in ^^N. 



employing R-matrix theory showed that the data could be best explained with a 100% 
sequential decay pattern. 



5.2.2. Two-proton emission from ^''Ne The decay of excited states of ^"^Ne by 2p 
emission was studied several times in recent years [HH ESI EH]. The first excited state of 
^^Ne is Ip bound but unbound with respect to 2p decay by 344 keV and may therefore 
decay by direct 2p emission to the ground state of ^^O (see Fig. [5]). Chromik et al. [8X1182] 
populated the first excited states of ^^Ne by Coulomb excitation and observed their decay 
by 7 de-excitation [81] or in a complete- kinematics experiment [82j . 

In the first experiment, the missing strength of the decay of the first excited state 
by 7 emission, as compared to theoretical calculations, was interpreted as a possible 
unobserved 2p branch. However, it was found that such an interpretation was in 
contradiction to barrier-penetration calculations. In fact, a 2p partial half-life shorter 
by a factor 1700 than predicted was needed to account for the missing decay strength. 
The second complete-kinematics experiment allowed determining all decay channels and 
a partial 2p half-life for the first excited state larger than 26 ps was deduced. No 2p 
emission from the first excited state was observed. The decay of the second excited 
state, also populated by Coulomb excitation, was observed to decay by sequential 2p 
emission to the ground state of ^^O yielding an isotropic emission pattern for the two 
protons. 

This experimental result was confirmed by Zerguerras et al. [83]. In their 
experiment, excited states of ^^Ne were populated by one- neutron stripping reactions 
from a ^^Ne beam at 36 MeV/nucleon. States up to 10 MeV excitation energy were 
observed to decay by 2p emission. The decay products were detected in the MUST 
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Figure 5. Decay scheme of the excited states in ^^Ne. The first excited state in -"^^Ne 
is bound with respect to one-proton emission but unbound with respect to two-proton 
emission to the ground state of ^^O. Hence, direct 2p decay is expected. For further 
details see the discussion in the text. 

detector array [S^j cind in the SPEG spectrometer [85j. These measurements allowed 
a reconstruction of the complete decay kinematics and therefore to determine the 
excitation energy of ^^Ne and the relative proton-proton angle in the center-of-mass 
frame. 
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Figure 6. Invariant mass spectrum for ^^O + 2p events ^83]. The peak at 18.5 MeV is 
interpreted as the decay of the second and third excited states of ^^Ne by two-proton 
emission, whereas the activity above 20 MeV arises from the decay of higher-lying 
states. 

Fig. [6] shows the invariant mass spectrum for ^^0-|-2p events. The angular 
distribution of the two protons was measured to be isotropic for the second and third 
excited states in ^'''Ne (Fig. Uh)- However, for higher lying states, a distinct angular 
correlation of the two protons could be observed for the first time (Fig. Uh)- This 
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result is rather surprising, as many intermediate states in the Ip daughter nucleus ^^F 
should be accessible and therefore a sequential decay mechanism via these states should 
dominate. 

The interpretation of this decay pattern is still controversial. Kanungo et al. [86] 
and Grigorenko et al. [87] suggested that some of the excited states of ^^Ne may have 
a pronounced 2p halo structure and a large overlap with the ^^O ground state leading 
to much larger spectroscopic factors for a direct 2p decay than for a sequential decay. 
Obviously, the validity of this assertion depends on whether the configuration mixing in 
^''Ne leaves the ^^O core intact. On the other hand, as shown in SMEC studies [33] the 
correlated 2p emission is always hindered if Ip decay channels are open. Moreover, the 
measured energies of the two emitted protons from ^^Ne show large differences which 
cannot be reconciled with a diproton scenario where roughly equal energies are expected 
for the two protons. 

Another explanation of the 2p emission pattern in ^''Ne could be the large 
deformation in higher lying states. Strong anisotropy of the Coulomb barrier could 
be a source of dynamical correlations between emitted protons even if this 2p decay is 
a sequence of two successive Ip emissions. Effects of this kind have been discussed in 
the de-excitation of heavy- ion residues (see e.g. Ref. |88J)- A similar explanation has 
been put forward also in connection with the anomalous decay properties of ^^Ag™ |27j . 
Obviously, for a consistent interpretation of the ^^O + 2p experiment, higher statistics 
data is mandatory. 
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Figure 7. Proton-proton angular distribution for 2p events with a ^^O recoil yielding 
a mass excess of less than 20.2 MeV (a) and above 20.2 MeV (b) (see Fig. E]). The 
solid and the dashed curves are simulations of a correlated two-proton emission via a 
^He resonance and a sequential emission pattern, respectively |83j . 

5.2.3. Two-proton emission from ^^Ne In an experiment performed at ORNL Oak 
Ridge, Gomez del Campo et al. [H9] used a radioactive ^^F beam on a (CH2)n target to 
produce excited states in ^^Ne. The experiment used a thick-target approach, i.e. only 
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the protons could be detected behind the target, whereas the beam-hke particles (e.g. 
the heavy ^^O recoils from the decay of excited states in ^^Ne) could not penetrate the 
target. Therefore, only protons with a relatively large energy loss depending on their 
emission point in the target could be observed. Using two different incident energies 
of the ^^F beam, these authors observed at higher incident energy the emission of two 
protons which they interpreted as a decay of the 6.15 MeV resonance state in ^^Ne. For 
the lower energy, where this state could not be populated, no coincident protons were 
observed. Both the relative proton-proton angle and the energy difference between the 
two protons was in agreement with simulations for a three-body decay of this state as 
well as for a diproton-type decay. This lack of sensitivity is to a large extent due to the 
limited statistics of the experiment, but also due to the limited angular coverage of the 
experimental setup. 

5.2.4- Two-proton emission from ^^Ag The decay of ^^Ag™ is unusual in several 
respects. It possesses two /^-decaying isomers and has the highest spin state ever 
observed which decays by (3 decay. In particular, the second isomeric level seems to 
decay by several different decay branches: /3-delayed 7 decay [901 El] , /3-delayed proton 
emission [92], direct proton decay [93], and direct 2p emission [27]. All these experiments 
were performed at the on-line separator of GSI [91] . 

The two protons were observed by means of silicon strip detectors which surrounded 
the implantation point of the on-line separator. This setup allowed for the measurement 
of the individual proton energies and the relative proton-proton angle. Although the 
statistics of the 2p experiment is rather low [27], the authors identified proton-proton 
correlations. They interpreted these data as due to simultaneous 2p emission from 
a strongly deformed ellipsoidal nucleus. In this case, the asymmetric Coulomb barrier 
favours the emission of two protons within narrow cones around the poles of the ellipsoid, 
either from the same pole or on opposite sides. 

One should mention, that the observed proton-proton correlations in such a scenario 
are largely unrelated to the proton-proton correlations inside a nucleus. The microscopic 
theoretical study of the ^^Ag*" decay is also impossible using the present-day computers, 
though only such studies could disentangle an internal structure of decaying states from 
dynamical effects of the anisotropic Coulomb barrier, providing more details about the 
2p decay pattern in the nucleus. Unfortunately, further experimental studies of the 
94^g"i (jgQay cannot be made without new technical developments to produce ^^Ag™, as 
the on-line separator of GSI has been dismantled a few years ago. 

6. Two-proton radioactivity 

Goldanskii [H [11] called it a pure case of 2p radioactivity, if Ip emission is not possible 
because the 2p emitter level is lower in energy than the level in the Ip daughter nucleus 
and, moreover, the corresponding levels are sufficiently narrow and thus do not overlap. 
The limitations of this picture will be discussed in chapter 8. Goldanskii also gave an 
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arbitrary half-life limit of 10~^^ s to distinguish between two-proton radioactivity 
and the decay of particle-unstable nuclei. None of the examples discussed above satisfy 
Goldanski's radioactivity criterion; they all have rather short half-lives (of the order of 
10^^^ s) and for the light ground-state emitters, the states which emit the first and the 
second proton overlap significantly. In this sense, 2p radioactivity was observed for the 
first time in "^^Fe experiments [T9| [20]. Later ^^Zn [21] and possibly ^'^Ni ^22j were shown 
to be other ground-state 2p emitters. 

Recent theoretical predictions [l5l[l7l[T6] agree that ^^Ti, "^^Cr, ^^Fe, and ^^'^^Ni are 
candidates to search for 2p radioactivity. Firstly, their predicted decay energy is close to 
or in the range of energies (1.0-1.5 MeV) where 2p radioactivity is most likely expected 
to occur, and secondly these isotopes are already or could soon be produced in projectile- 
fragmentation type experiments. This production technique allows to unambiguously 
identify the isotopes before their decay and thus measure the decay half-life and the 
decay energy, the necessary ingredients to identify this new radioactivity. The most 
promising candidates were "^^Fe and "^^Ni with decay energies of 1.1-1.4 MeV. 

6.1. Identification of the 2p candidates and first negative results 

The first of these isotopes, ^^Ti, was observed at the Grand Accelerateur National Ions 
Lourds (GANIL), Caen [96] in 1990. This experiment was also the first attempt to 
observe 2p radioactivity of these medium-mass nuclei after projectile fragmentation. 
However, no indication for 2p radioactivity could be found. Instead, a /5-decay half-life 
of 28 l^ms was measured and /^-delayed protons were observed. Another attempt to 
search for 2p radioactivity from ^^Ti was made at GANIL in 1991 with an improved 
setup [971 [98]. However, again no 2p events could be observed. 

The nuclei *^^Cr, *^^Fe, and ^^Ni were identified for the first time a few years later 
at GSI ^9]. In this experiment, no spectroscopic studies could be performed and 
the observation of these isotopes was the only outcome. Nonetheless, this experiment 
provided the first observation of one of the most promising 2p radioactivity candidates 
and opened the door for the search of ^^Ni, another candidate. This doubly-magic 
nucleus was then discovered in 1999 [100] in a GANIL experiment using the SISSI/LISE3 
device [101] . In the same experiment, the observation of ^^Fe was confirmed with about 
a factor of 10 higher statistics. 

The GANIL experiment was also set up to perform spectroscopic studies of the 
very exotic nuclei produced [102] . The decay trigger in this experiment allowed mainly 
/9-delayed events to be registered. /5-delayed decay events from e.g. ^^Fe decay had 
a detection efficiency of about 30%, as they were triggered by two detectors adjacent 
to the implantation detector. 2p events without (3 particle had a very low trigger rate 
(a few percent) due to a heavy-ion trigger for the implantation device. This situation 
allowed for the observation of the decay of "^^Fe, however, basically of only its /3-delayed 
branch (see Sect. 6.2). In any case, no claim about the observation of 2p radioactivity 
was made. 
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Nevertheless, the experiment allowed one to test two other possible, although less 
promising, 2p candidates: ^^Cr and ^^Ni. As mentioned above, these two nuclei were 
predicted to be possible 2p emitters, however, with decay energies most likely too small 
to exhibit this decay. The decay energy spectrum obtained during this experiment for 
^^Cr is shown in Fig. [8l The prominent peak at 1.9 MeV was excluded to be of 2p 
nature, as such a high energy would lead to an extremely short barrier penetration time 
for the two protons (see Fig. [HI right-hand side), in contradiction with the ^^Cr half-life 
of 13.4t2!4ms measured in the same experiment. 




Figure 8. Left-hand side: Charged-particle decay-energy spectrum obtained after 
^^Cr implantation. The 1.9 MeV peak is interpreted as a /3-delayed decay |102j . 
Right-hand side: Barrier penetration half-hfe for a di-proton in ''^Cr as a function 
of the available decay energy [102j . The red line represents the 1.9 MeV peak from the 
decay-energy spectrum yielding a barrier penetration half-life of 10~^^s if it were of 2p 
origin. 



In a similar way, ^^Ni could be excluded from having a significant direct 2p branch. 
The last 2p emitter candidate tested in this experiment was ^^Ti. In agreement with 
previous experiments, this nucleus turned out to decay predominantly by (3 decay. 
Several /3-delayed proton groups could be identified and the 2p emission Q value (Q2p = 
670±100 keV) could be determined, indicating a rather high barrier-penetration half-life 
and therefore a rather low probability for 2p radioactivity. 



6.2. The discovery of two-proton radioactivity 

The discovery of ground-state 2p radioactivity had to await the year 2002, when two 
experiments, one performed in 2000 at the SISSI/LISE3 facility of GANIL [19] and one 
performed at the FRS of GSI in 2001 [20j, were published. Both experiments used the 
fragmentation of a ^^Ni primary beam. The fragments of interest were implanted in 
a silicon detector telescope (see Fig. [9]), where correlations in time and space between 
implantation events and decay events could be performed. This technique allowed to 
unambiguously correlate "^^Fe implantation events with their decays. 
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Figure 9. Left-hand side: Detection setup used in the GANIL experiment to 
identify two-proton radioactivity. The setup consisted of a sihcon telescope, where 
the implantation device was a double-sided silicon strip detector (DSSSD). Several 
time-of- flight, energy-loss and residual energy measurements allowed an unambiguous 
identification of the projectile fragments. The radioactive decay of the isotopes of 
interest was detected in the DSSSD. The last detector, a Si(Li) detector, was used 
to detect /? particles [19]. Right-hand side: Setup used in the GSI experiment. A 
silicon telescope of seven elements was surrounded by a high-efficiency Nal array. 
Implantation and decay events were correlated in different detectors depending of the 
implantation depth of the ^^Fe ions. The Nal device served to search for 7 radiation 
from positron annihilation |20j . 



Fig. [To] shows the decay energy spectra from the two experiments. Both find a 
decay energy of about 1.1 MeV with a decay half-hfe in the 3-5 ms range. An additional 
feature of both experiments was the capabihty to detect additional radiation like /3 
particle in the GANIL experiment, or 511 keV annihilation radiation and other 7 rays 
in the GSI experiment (see Fig. [9]). This capability allowed for a rejection with a high 
probability of a possible /3-delayed origin of the decay events with an energy release 
of about 1.1 MeV. Other pieces of evidence from the GANIL experiment were the 
smaller width of the 1.14 MeV peak as compared to /3-delayed proton emission peaks in 
neighbouring nuclei and the observation of daughter decays consistent with the decay of 
^^Cr, the 2p daughter of ^^Fe. In addition, the experimentally observed decay energy is 
in nice agreement with theoretical predictions flE l [T7 1 [T6]. These experimental results 
led the authors to conclude the discovery of ground-state 2p radioactivity in the decay 
of ^^Fe. 

Of particular interest was also the observation of decay events from ^^Fe with an 
energy release different from the 1.1 MeV of the main peak. Some of these additional 
events could be identified as being in coincidence with (3 particles and were therefore 
attributed to /^-delayed decay branches. From these first results, a 2p radioactivity 
branching ratio of 70-80% was determined for ^^Fe. The partial half-lives of ^^Fe for 
2p emission and for j3 decay are therefore of the same order and allowed the two decay 
modes to compete. 
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These results on ^^Fe have recently been confirmed in a new GANIL experiment |22j . 
Several experimental improvements allowed a better energy calibration and a smaller 
data acquisition dead time. The decay energy was determined to be 1.154(16) MeV, 
the half-life to be 1.6l[5:| ms, and the 2p branching ratio to be 57(10)%. Together 
with the results from the two earlier experiments [191 ED], this experiment allowed 
to determine average values for the decay energy (Q2p= 1.151(15) MeV), the half-life 
(Ti/2= 1.75lo28ms), the 2p branching ratio (BR = 59(7)%), and the partial 2p half-life 
{Tl%= 3.0t°|ms). 

In 2007, a new experiment performed at the A1900 separator of Michigan State 
University allowed the observation of 2p events in a TPC (see below). The newly 
determined half-life and the 2p branching ratio are included in Table [1] to determine 
average values. These average values will be used in the comparison with theoretical 
model of 2p radioactivity. 
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Figure 10. Left-hand side: ^^Fe decay-energy spectrum from tire GANIL experiment 
exliibiting a peak at (1.14±0.04) MeV [T^. Right-hand side: "^^Fe decay-energy 
spectrum from the GSI experiment showing four events at 1-1(1) MeV [20]. Both 
spectra also show events at higher energies. 

Another outcome from the GANIL experiments [T9l [22] on ^^Fe was the 
determination of the half-life of the 2p daughter nucleus. From the observation of the 
second decay after ^^Fe implantation, the half-life of the daughter decay, in coincidence 
with the 1.151 MeV peak, could be determined. In Fig. [11], this half-life is compared to 
the half-lives of all possible "^^Fe daughter nuclei. Only the half-life of "^^Cr is in agreement 
with the experimentally observed daughter half-life, thus giving an independent proof 
for the observation of 2p radioactivity. 
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Figure 11. The experimentally observed half-life for the decay of the daughter of ^'^Fe 
is compared to the half-lives of all possible daughters [73]. Only the half-life of ^^Cr is 
in agreement with the experimental value. 



6.3. Two-proton radioactivity of Zn 

The nucleus ^^Zn was early identified as a possible 2p emitter [101 [H]. More recent 
Q2p estimates yielded decay energies of 1.79(12) MeV [T6], 1.87(24) MeV [I03], and 
1.33(14) MeV |104j . Due to rather large theoretical error bars and the uncertainty 
concerning the decay mechanism, it was rather unclear whether ^^Zn would exist at all, 
i.e. whether its half-life would be long enough to be observed in projectile-fragmentation 
type experiments. For this purpose, a life time of the order of a few hundred nano- 
seconds was mandatory. 

The search for ^^Zn began with the observation of two other proton-rich zinc 
isotopes, ^^Zn and ^^Zn |105] . They were observed for the first time in experiments 
at the LISE3 separator |101] via 2p pick-up reactions with a ^^Ni primary beam. The 
identification of these two isotopes allowed also for an extrapolation of the production 
cross section of ^^Zn. 

The nucleus ^^Zn was synthesized and observed for the first time in 2004 in an 
experiment at the SISSI/LISE3 facility of GANIL Caen [21]. Similarly to the experiment 
with "^^Fe, ^'^Zn was produced with a primary ^^Ni beam at 75 MeV/nucleon impinging 
on a natural nickel target. The fragments of interest were selected and separated with 
the LISE3 separator and implanted in a silicon telescope with a double-sided silicon 
strip detector being the central device. The ^^Zn nuclei were identified by means of 
time-of-flight, energy- loss and residual-energy measurements. Fig. [12] shows the two- 
dimensional identification plot for the ^^Zn experiment. All in all, eight ^^Zn nuclei 
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could be unambiguously identified. 
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Figure 12. Two-dimensional identification plot for ^"^Zn and neighbouring nuclei. The 
energy loss in the first silicon detector of the detection setup is plotted as a function 
of the time-of-flight through the separator. Additional parameters are used to purify 
the spectrum [21]. 

The setup allowed correlating in time and space these implantation events with 
subsequent decays. Seven of the eight implantation events are followed by a decay 
with an energy release of 1.48(2) MeV (see Fig. [T3l) . None of these decay events is in 
coincidence with a f3 particle. As for ^^Fe, the 1.48 MeV peak is much narrower than 
/3-delayed proton lines in neighbouring nuclei. Although less significant as in the case 
of ^^Fe, daughter decays in agreement with the decay of ^^Ni, the 2p daughter of ^^Zn, 
could be observed. Finally, the experimental half-life (Ti/2= 3.2lo]8^^s) is in agreement 
with what is expected for a 2p decay with an energy release of 1.48 MeV. From these 
experimental observations, the authors concluded the observation of ground-state 2p 
radioactivity from ^^Zn. The 2p branching ratio was determined to be 871^7%, which 
yielded a partial half-life for 2p decay of 3.7li;o ^s. These results will be compared to 
theoretical predictions in Sect. 9.3. 

6.4- The decay o/^^Ni 

The discovery of ^^Ni and the observation of its decay was considered important for 
several reasons. Firstly, it is a doubly-magic nucleus which was predicted to be particle 
stable or quasi-stable only due to strong shell effects. Secondly, ^^Ni is one of the most 
exotic of all doubly-magic nuclei within experimental reach. Its observation and the 
measurement of its properties, like e.g. the excitation energy of the first 2"*" state, were 



Two-proton radioactivity 



22 



(0 



o 

O 




1000 1500 



2000 2500 3000 3500 4000 

Energy (keV) 



Figure 13. Charged-particle spectrum from the decay of ^^Zn. The peak at 
1.48(2) MeV is due to two-proton emission from the ground state of ^^Zn. The event 
at 4.19 MeV is in coincidence with a /3 particle in the adjacent detector [5T]. 



expected to teach us a lot about the persistence of shell structure far from the valley of 
stabihty. 

The nucleus ^^Ni was rather early identified as a possible 2p emitter [1]. This 
prediction was confirmed and refined by new mass and Q2p calculations, which yielded 
1.36(13) MeV [15], 1.14(21) MeV [Hj, 1.35(6) MeV [16], and 1.29(33) MeV [103]. When 
compared to ^^Fe and ^"^Zn, it is evident that this Q2p value is again in the right range 
to permit "^^Ni to have a significant 2p radioactivity decay branch. 

The isotope ^^Ni was first observed in 1999 at GANIL Caen [100] . where four events 
could be unambiguously identified. This result was confirmed in 2004, when again four 
events of "^^Ni could be registered in a new SISSI/LISE3 experiment at GANIL |22j . 
These data were obtained in the same run as the second ^^Fe data set from GANIL, 
because the acceptance of the SISSI/LISE3 device allows the measurement of different 
neighbouring nuclei at the same time. Fig. [Til shows the fragment identification plot. 

The four ^^Ni implantation events could be correlated in time and space with 
subsequent decay events. From the four decay events, only one event has all 
characteristics of a 2p emission: a decay energy of 1.35(2) MeV in the expected energy 
range and no coincident f3 particle. The decay takes place after 1.66 ms and the daughter 
decay is also fast as expected for "^^Fe, the 2p daughter of "^^Ni. These four decay events, 
one of possible 2p origin and three from /3-delayed decays, were also used to determine 
the half-life of ^^Ni (Ti/2= 2. lloiyins). Although this single possible 2p decay event is 
not considered sufficient to claim the observation of 2p radioactivity for ^^Ni, these data 
were nonetheless used to determine the 2p branching ratio to be 251^9% ^^e partial 
2p half-life to be Tf',^=8At^f^^ms [22]. 
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Figure 14. Two-dimensional identification plot for ^^Ni and neighbouring nuclei. The 
energy loss in the first silicon detector is plotted as a function of the flight time through 
the separator. Additional parameters are used to purify the spectrum [22] , 



6.5. Direct observation of two protons in ground-state decay of^^Fe 

From all 2p emitters observed up to now, the ground-state emitter ^^Fe is the best 
candidate to search for a clear signal of correlated emission and to study in detail the 
decay mechanism. 

As just mentioned, 2p emission could only be governed by phase space or there could 
be more or less strong angular and energy correlations. In order to study these questions, 
new experimental setups had to be imagined, the basic limitation of the silicon-detector 
setups being that the 2p emitters are deeply implanted in a silicon detector, which then 
does not allow the detection of individual protons, but rather of the total decay energy 
and the half-life only. 

The solution to this problem is the use of gas detectors capable to visualize the 
traces of the emitted protons. Therefore, time projection chambers (TPC) were set 
up in Bordeaux |106j and Warsaw [107] . The Bordeaux detector is based on gas 
electron multiplier technology |108j and a double-sided microgroove detector |109j with 
an application-specific integrated-circuit read-out, whereas the Warsaw detector uses 
an optical read-out system by means of a digital camera and a photomultiplier tube 
(PMT) which registers the time sequence. Both TPC detectors have recently taken first 
data [721 IllOj (see Figs. [15], [16]). Due to the highest production rates, "^^Fe was the first 
2p emitter to be studied with these devices. 

Due to the higher statistics, the data taken at Michigan State University with the 
Warsaw TPC allowed to compare the experimental data to predicitons of the Si'^M 
model of Grigorenko et al. model [lllj . for which nice agreement for the angular and 
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Figure 15. A single event of two-proton decay of ^^Fe from a recent experiment 
performed with the Bordeaux TPC at the LISE3 facihty of GANIL P72]. The top row 
shows the "^^Fe implantation event, where the ion enters the chamber parallel to the X 
direction and stops in the center of the chamber. The decay (bottom) takes place at 
the point where the ion stopped. 

energy distribution was reached [110] indicating e.g. an equal sharing of the energy 
between the two protons. 

7. Conclusion for the experimental part 

The observation of ground-state 2p radioactivity for ^^Fe, ^^Zn and possibly for '^^Ni 
allowed this type of radioactivity to be established as a new nuclear decay mode. 
Although only based on theoretical mass predictions, these decays are believed to be 
simultaneous emissions of the two protons. A sequential decay could have a sizable 
contribution only if the Qip values were off by a few hundred keV what is rather 
unlikely. Nonetheless, the assumption that the decay is indeed simultaneous still awaits 
experimental verification though, as will be discussed in Sect. 9, a symmetric energy 
distribution of individual protons emitted in the 2p decay of ^^Fe |110j is a strong 
indication for a predominance of the simultaneous 2p emission. Table [1] summarizes the 
experimental information available for ^^Fe, ^^Ni, and ^^Zn. 

To get a coherent picture of this new decay mode, more 2p emitters have to be 
studied experimentally. However, as these nuclei are, by definition, situated beyond the 
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Figure 16. An example of a two-proton decay event of ^^Fe recorded by the Warsaw 
TPC detector. Left: an image taken by the CCD camera in a 25 ms exposure. The faint 
track of a ''^Fe ion entering the chamber from the left is seen. The two bright, short 
tracks are low-energy protons emitted approximately 2.5 ms after the implantation. A 
weaker and longer track, going from the decay vertex to the lower right, represents a 
high-energy proton emitted about 8 ms after the implantation which is due the decay of 
''^Cr, the 2p decay daughter of ''^Fe. Right: the time profile of the total light intensity, 
measured by the PMT, provides the time sequence of the decay events. In the inset, 
a magnified part is plotted showing in detail the 2p emission. Figure courtesy of M. 
Pfiitzner. 
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Table 1. Decay energies, half-lives, branching ratios, and partial half-lives as 
determined in three independent experiments for the decay of '*^Fe [I9 l l20 1 f22 llll0j as 
well as for ^SNi ^22^, and ^^Zn [21] • 



proton drip line, their production is rather difficult and the rates are very low. Rather 
high rates can be achieved today at the SISSI/LISE3 facility of GANIL with about 
10-15 ^^Fe per day, about 2 ^^Zn per day, and roughly 1 ^^Ni per day under optimal 
conditions. It was recently shown that at the coupled-cyclotron facility and the A1900 
separator of the NSCL even higher rates can be achieved. Other candidates like ^^Ge, 
^^Se or ^^Kr might be reached in the near future at the new RIKEN Radioactive Ion 
Beam Facility (RIBF) in Tokyo, Japan. 
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8. Theoretical description of two-proton radioactivity 

The nuclear many-body Hamiltonian is supposed to describe all nuclei that can exist 
and not merely one nucleus of a given number of protons and neutrons. In this sense, a 
nucleus is never a closed, isolated quantum system but communicates with other nuclei 
through virtual excitations, decays and captures. The communication is broken and the 
nucleus becomes artificially closed if the subspace of continuum states is excluded in 
the network of coupled systems. Obviously, the closed quantum system idealization of 
a real many-body system has very different features from those observed [32| 1112] . in 
particular in the neighborhood of each reaction threshold |113[ 1114] . 

Most of the nuclear states are embedded in the continuum of decay channels 
due to which they get a finite lifetime. The initial impulse for a mathematical 
formulation of the shell model (SM) for open quantum systems goes back to Feshbach 
and Fano |115l I116[ 1117] who introduced the two subspaces of Hilbert space: (i) 
the subspace of discrete states, and (ii) the subspace of real-energy scattering states, 
respectively. Feshbach achieved a unified description of direct processes at the short- 
time scale and of compound nucleus processes at the long-time scale |115l 1116] . A 
unified description of both nuclear structure and nuclear reactions in a single theoretical 
framework became possible only in the continuum SM [118\ 11191 1120^ 11211 1122^ 1123] . 

First realistic studies using the continuum SM have been presented only at the 
turn of the century in the framework of SMEC [301 EI] (see also Ref. [32] for a 
recent review). This model is based on the completeness of the single-particle basis 
consisting of bound orbits and a real-energy scattering continuum. Single-particle 
resonances are regularized, i.e. they are included in a discrete part of the spectrum 
after removing scattering tails which are incorporated in the embedding continuum. 
Subsequently, the new set of discrete states and the non-resonant scattering states are 
reorthogonalized. Configuration mixing in the valence space (the internal mixing) is 
given by an effective SM interaction. Coupling of discrete states to the continuum of 
decay channels is calculated microscopically and leads to the so-called external mixing 
of SM wave functions. The decay channels are obtained in the S'-matrix (scattering- 
matrix) formalism (see Ref. [32] and references quoted therein). 

Another approach referring to the continuum SM has been proposed recently in 
Refs. [12411125] . Internal mixing in this approach is given by an effective SM Hamiltonian 
but an external mixing is neglected. This scheme has been applied to the description of 
observables in the one- and two-particle continuum |125] . 

The genuine three-body decay models for 2p radioactivity have been proposed in 
Refs. |126l 11271 11281 1129] using the hyperspherical harmonics method. These models 
combine the cluster model description of the decaying system with a treatment of 
three-body asymptotic states of charged particles. Rigorous treatment of a three- 
body kinematics in these models allows to study both angular and energy /momentum 
correlations between decaying fragments. 

A different approach for describing many-body open quantum systems is provided 
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by a generalization of SM in the complex- momentum plane (the GSM) |130l [551 1131j . 
This general formalism with no restriction on the number of particles in the scattering 
continuum has not yet been applied for a description of experimentally observed 2p 
decays. 

8.1. Shell Model Embedded in the Continuum with one particle in the scattering 
continuum 

A general theory of the decaying many-body system is provided by the open quantum 
system formulation of the SM. Before we turn to the description of 2p decay, let us 
introduce the basic features of SMEC in the simpler case of the one-particle continuum. 
In SMEC, the function space consists of two sets: the space of £^-functions, as in the 
nuclear SM, and the space of scattering states (decay channels) embedding SM states. 
The decay is a result of the coupling between discrete states and decay channels. These 
two sets of wave functions are found by solving the Schrodinger equation for discrete 
states: 

HsMm=Ef''^\<l>,) , (1) 
where Hsm is the SM Hamiltonian, and for scattering states: 

Y^iE - H^,)^i'^^^ = , (2) 

c' 

where ifcc = Ho + Vcc in (El) is the standard coupled-channel Hamiltonian. Different 
channels are defined by coupling the motion of an unbound particle relative to the 
residual nucleus with (A — l) particles in a discrete SM state |$/"^). are scattering 
states projected on the channel c, and the channel numbers c are defined by the quantum 
numbers of of the states j of the residual nucleus and those of the unbound particle which 
are coupled to the total quantum numbers J'^ and T of the system. 

Using these two function sets: Q = V = {C,e}, one defines projection 

operators: 

Q = j:\^t){^t\ ; P= dE\^E){U (3) 

i=l -^o 

satisfying Q\^e) = 0, P\^f) = 0, and projected Hamiltonians: QHQ = Hqq and 
PHP = Hpp. The closed system Hamiltonian Hqq is the nuclear SM Hamiltonian 
Hsm, whereas Hpp is given by Hcc- To construct a complete solution in Q + V = Id'- 

\^e) = Q\'^e)+P\^e) , 

one has to find a wave function connecting Q and V subspaces. This wave function cjf^'* 
is obtained as a solution of non-homogeneous coupled-channel equations: 

\u;t\E)) = G^;\e)\w,) , (4) 

where \wi) = Hpql^i) is the source term, 

G^p\e) = \im[E + ze- Hpp]-^ 
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is the Green's function in V subspace, E is the total energy of the nucleus [A], and 
HpQ = PHQ. Using the three function sets: {^e}, and ujI~^\e), one can find 

the solution in the total function space: 

I^^e) = \^e) + E(l*f ) + \ujt\E)mt\{E - nQQ{E))-'\^t){^^\HQp\^j,) . (5) 

i,k 

1-Lqq{E) is the energy-dependent effective Hamiltonian in the function space of discrete 
states: 

Hqq{E) = Hqq + HqpG^^\e)Hpq (6) 

which takes into account couplings to decay channels. In that sense, TCqq is the open 
quantum system Hamiltonian in Q subspace. 

The second term on the right-hand side of Eq. ([6]) provides an imaginary part 
of energy eigenvalues for states above the particle emission threshold. The real part 
of this term contributes to the real part of energy eigenvalues both below and above 
the threshold. For that reason, it would be more appropriate in phenomenological 
applications to fit matrix elements of Tigg and not of Hqq to nuclear spectra in a given 
mass region. This ambition task has not yet been accomplished. 

Energies Ei and widths Fj of the resonance states derive from the solutions of 
fixed-point equations: 

E, = E,{E = E,) ; Ti = ti{E = Ei) (7) 

where functions Ei{E) and Ti{E) follow from (complex) eigenvalues of H-qq. 

The eigenvalues of TCqq have a physical meaning if the total function space is divided 
into the subspaces of the system (the subspace Q) and the environment (the subspace 
V) as follows: the system contains all resonance-like phenomena while the environment 
describes the smooth part in the energy region considered. This operational definition 
of subspaces implies that in SMEC one has to construct carefully the single-particle 
basis. The narrow single-particle resonances in the continuum are included into the set 
of discrete states either as resonance anamneses [132\ or as quasi-bound single-particle 
states [133] . The new definition of discrete single-particle states implies a redefinition 
of the continuous spectrum, as described in Ref . [32] . 

8.2. Theory of two-proton radioactivity 

The description of Borromean (halo) systems and two-nucleon decays involves couplings 
to the two-nucleon decay channels. In the 2p decay, since two protons do not form a 
bound two-particle system, one neither can uniquely specify the final state of two protons 
nor the transition path leading to this state. Various transition paths may interfere 
making distinction between different decay mechanisms difficult if not impossible. In 
this respect, the 2p decay is qualitatively different from any standard two-body decays 
such as the proton-, neutron- or a-particle decays. 

A generalization of the SMEC formalism for larger number of particles in the 
scattering continuum is conceptually similar to the one-particle case discussed above 
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(cf Sect. 18. ip . For the description of two-particle continuum, one has to introduce a 
subspace T with two particles in the continuum, and T the corresponding projection 
operator. The completeness relation in total function space is then: Q + V + T = 1. 
Similarly as in Sect. 18. one can decompose the Hamiltonian H into components acting 
in different subspaces and the coupling terms between those different subspaces. The Q 
subspace effective Hamiltonian Hqq{E) can be written as [31]: 

'^Qq{E) = Hqq + HQpGj.{E)HpQ 

+ [Hqt + HqpG%{E)Hpt\G^{E)[Htq + HtpG%{E)Hpq\ (8) 

The first term on the right-hand side is the closed quantum system Hamiltonian, the 
second term describes a coupling with the one-particle continuum, and the third term 
corresponds to all possible couplings with the two-particle continuum. G^{E) is the 
Green function in T 

G^{E) = \im[E + ie- Htt - HtpG%{E)Hpt\~^ 

modified by couplings to V. Similarly as in the case of one-particle continuum, one 
defines the source term 

\wi) = [Htq + HtpG^^\e)Hpq]\^,) (9) 

and the continuation \uj^) of SM wave function into the two-particle continuum: 

Icot) = G^^\e)[Htq + HtpG^^\e)Hpq]\^,) = G^jt\E)\w,) . (10) 

Matrix elements of the operator responsible for Q — T couplings (the third term on the 
right-hand side of Eq. ([8])) can be expressed as the overlap between the source term 
and the function Uj' (see Eqs. ([9]) and ( |TOl) ). 

Effective Hamiltonian T-Cqq{E) can be written in various equivalent forms. For our 
purpose, it is convenient to extract the direct coupling term between Q and T subspaces: 

'^Qq{E) = Hqq + HqtG''j^\e)Htq 
+ [Hqp + HQTG^rj^\E)HTp]G^p\E)[HpQ + HpTG^rj^\E)HTQ] . (11) 
Gp^^ (E) in (fTTj) stands for P-subspace Green function modified by couplings to T 

&p\e) = \im[E + ie- Hpp - HptG^^\e)Htp\-^ , 

e— >0 

where 

G^^\e) = lim [E + ie- HttT^ 

is the Green function in T. 

Hqq{E) (cf Eqs. ([8]) and (ITT]) ) describes all emission processes involving one and 
two nucleons (two protons). Observed 2p partial decay width is a sum of contributions 
of different interfering processes which cannot be disentangled experimentally. In 
theoretical studies, one can always isolate a given processes by switching-off all others 
and calculate the contribution of such an isolated process to the 2p partial width. 
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However, relative weights of different isolated processes do not include interference 
effects and, hence, the resulting picture of the 2p decay is somewhat distorted. In 
this sense, and contrary to the two-body decays, the interpretation of experimental 2p 
decay data is unavoidably model dependent. 

In the observed 2p decays, certain emission scenarios may be less probable, so it 
is instructive and certainly less cumbersome to consider limiting cases of the 2p decay. 
For example, if V ^ T couplings are weak then the dominant process is a direct 2p 
emission described by 

n'^S\E) = Hqq + HqpG^^\e)Hqp + HqtG^^\e)Htq , (12) 

where the second term on the right-hand side is responsible for couplings to one-nucleon 
decay channels. On the contrary, if an intermediate system [A — 1] plays an important 
role, then Q ^ T couplings can be neglected and the effective Hamiltonian reads: 

'^qq\e) = Hqq + Hqp&p\e)Hpq . (13) 

This operator describes a standard sequential 2p emission via an intermediate resonance 
in the [A — 1] nucleus. 

The sequential 2p emission process may also occur via continuum states which are 
correlated by either weakly bound or by weakly unbound states of [A — 1] nucleus. This 
process which goes through the tail of an energetically inaccessible state, is called the 
virtual sequential two-body decay |134] . For its study, it is important to separate Ip- 
and 2p-decays in the effective Hamiltonian: 

'^Qq\E) = Hqq + HqpG^p\E)Hpq 

+ Hqp&^^HptG^^\e)HtpG^;\e)Hpq . (14) 

The second term in Eq. (IT^ describes the Ip decay, whereas the third term describes 
the genuine sequential decay. 

8.2.1. Virtual sequential two-proton emission The sequential 2p decay via the 
correlated scattering continuum of [A — 1] nucleus (the virtual sequential emission 
process) is an important component of the 2p decay mechanism. An effective 
Hamiltonian for this process ( cf Eq. ( fT4l) ) describes interference effects between the 
Ip decay and the sequential 2p decay, both for closed (virtual Q-V excitations) and 
opened (Ip decays) Ip decay channels. Diagonalizing T-Cqq\E), one obtains energies 
of states in [A] nucleus as well as their widths associated with the emission of one 
and two protons. In general, these two decay modes cannot be separated one from 
another. However, since couplings responsible for the 2p decay are smaller than those 
associated with the Ip decay, one may proceed in two steps. In the first step, one 
neglects the term HqpG^p^ HpTGiit\E)HTpG'-p\E)HpQ (cf Eq. ([14])) responsible for 
the 2p decay, and diagonalizes the remaining operator in a SM basis {1$'^)}. This 

~ A 

provides new basis states {|$ )} which are linear combinations of SM states. In the 
second step, using those mixed SM states one calculates the sequential 2p emission width 
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of a decaying state i.e. one calculates an imaginary part of the matrix element 

{^f\HQpGpHpTG^^\E)HTpG'^+\E)HpQ\^f) . 

Let us assume that different steps of the sequential process are uncorrelated, i.e. 
the first emitted proton is a spectator of the second emission. This implies: 

Hpp — > Hq/q' + p/iop 
Hj"p — > Hp' pi + p/ioP 

where primed quantities refer to {A — l)-nucleon space. In the Q' subspace, {A — 1) 
nucleons are in quasi-bound single-particle orbits, whereas in the V' subspace, one 
nucleon is in the continuum and {A — 2) nucleons are in quasi-bound single-particle 
orbits, ho is a one-body potential describing an average effect of [A — l)-particles on the 
emitted proton and p denotes a projector on one-particle continuum states. With this 
identification, HpT becomes a coupling between newly defined Q' and V subspaces. 

The virtual sequential 2p emission appears always if Q2p > 0. In that sense, the 
virtual sequential process is a component of any quantum-mechanical three-body decay. 
In the approximate scheme employed in SMEC, the virtual sequential 2p decay process 
competes with the diproton mode even for closed Ip emission channels. 



8.2.2. Diproton emission The effective Hamiltonian (fT2l) describes both the three- 
body direct 2p decay and its limit of two sequential two-body decays. In this section, 
we shall discuss the latter two-step scenario [29]. In the first step, protons are emitted 
as a '^He' -cluster which, in the second step, disintegrates outside of a nuclear potential 
due to the final state interaction |135[ 1136] . The final state proton-proton interaction is 
taken into account by the density p{U) of proton states corresponding to the intrinsic 
energy U in the proton-proton system pL37J. The total system [A — 2] (g) [2p] is described 
as the two-body system in a mean-field Uq: 



Htt = f^"^^ 



2/i 



f("') , (15) 



where jg ^j^g intrinsic Hamiltonian of the daughter nucleus [A — 2], and H^'^'-'> 

is the intrinsic Hamiltonian of '^He'-cluster. The cluster is described as a particle of 
mass M = 2Mp {Mp denotes the proton mass) and charge Z = 2. P^/2/i is the intrinsic 
kinetic energy of the system [A — 2] ® [2p], and fi is the reduced mass of the system. 
T*^'^') is the projection operator on the subspace of cluster states in the continuum of the 
potential P^/2/i -|- Uq- The 2p decay width is given by the matrix element: 

Tpp] = - 2/m {{^t'^\HQTG^^'-\E)HrQ\^t'^)) (16) 

where l^f"*"*) is the intrinsic state corresponding to a mixed state in the parent 
nucleus, {wf) = i^Tgl^l*"^*'*) is the source term, and \ujJ''^~^^) = G^\E)HTQ\4>i'"'^^) is a 
continuation of in T. The Coulomb interaction is included as an average field 

and does not enter in Htq- The matrix element (fT6j) is the integral over an intrinsic 
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energy U of the overlap of a channel-projected source term and a channel-projected 
continuation in T of an intrinsic state Ic^f""*^), weighted by the proton state density 
p{U) [33l[3l]. Working with intrinsic states allows to account for the recoil correction. 

The decay energy is divided into an intrinsic energy of the '^He'-cluster and its 
center of mass motion energy according to the proton-state density [33l [3l]. In this 
way, intrinsic degrees of freedom of the two protons are reduced and, therefore, the 
three-body decay is reduced to the two sequential two-body decays. 

K-matrix model of Brown and Barker This model extends the R-matrix approach for 
SM description of a three-body decay |137[ [29] . Similarly as in the SM approach for 
two-body decays, one may assume that the three-body decay is independent of the small- 
distance many-body dynamics. This opens for a definition of three-body spectroscopic 
factors which can be computed in the SM framework [29j in an analogous way as the 
preformation factors used to describe the a-particle decay. 

The R-matrix model of a diproton decay |137l [29] includes the s-wave proton- 
proton interaction as an intermediate state. This final state interaction is crucial for a 
quantitative description of experimental data, changing the diproton lifetimes by several 
order of magnitudes. ^He preformation factor is taken into account by ^He spectroscopic 
factors S[2p]. To calculate them, the SM wave function is projected onto the Os internal 
relative wave function for a '^He '-cluster |29j : 

S[2r>]={-^)'G'CiA,Z). 

In /p shell model sace, = 5/16, A = 6. C{A,Z) = | (^(A - 2.Z - 2)|^c|*(A, Z)) p 
is the overlap for the ^He wave function \E'c with L = 0, S = 0, and T = 1 in the SU{3) 
basis. 

The extended R-matrix model [29] neglects both the small- distance dynamics and 
the virtual scattering to continuum states. The latter process changes the ^He emission 
probability by modifying the wave function of a decaying state and changing the ^He 
preformation factor [33j. This change is expected to be particularly strong for states 
close to the 2p emission threshold(s), e.g. the ground states of ^^Fe, ^'^Ni and ^^Zn. 
Reasons for the continuum-induced change of ^He preformation probability or one- 
nucleon spectroscopic factors [1381 1139[ 1140 ] are the same as for a dominance of cluster 
states close to their respective decay thresholds [141J : both phenomena are consequences 
of the 'alignment' of near-threshold states with the channel state [32 | 11421 11431 1144[ 1113] . 

It is not easy to assess quantitative consequences of those simplifications in the 
extended R-matrix model. In the following sections, we shall try to address this problem 
by (i) comparing SMEC results for diproton decays with those of an extended R-matrix 
approach, and (ii) by analyzing SMEC results for virtual sequential 2p emission with 
and without external mixing. 

8.2.3. Direct 2p emission with three-body asymptotics In general, structures in the 
continuum are more difficult to treat than bound-state problems. Firstly, the three- 
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body problem for unbound systems is much less established, although studied for short 
and long-range interactions |145l 11461 135]. Secondly, the Coulomb problem with three- 
body asymptotics is still considered unsolved |147[ 11481 1149] even though an important 
progress has been reported recently in the description of proton-deutron scattering and 
of three-nucleon electromagnetic reactions involving ^He |150j . An infinite range of the 
Coulomb interaction in channel-coupling potentials does not allow to decouple equations 
at infinity. Consequently, an asymptotic behaviour of ujj'^J{p) cannot be defined without 
a certain approximation. One approximate way to proceed is to neglect off-diagonal 
channel-coupling potentials above a certain value of the hyperradius and to define an 
effective Sommerfeld parameter from diagonal coupled-channel potentials |128] . This 
technique has been applied in the three-cluster model of the 2p emission |15H 11281 1152j . 

In order to calculate the direct 2p emission probability, one has to evaluate the 
matrix element ( |T6i) without a simplifying assumption of two sequential two-body 
decays. Even though this three-body problem has been conveniently formulated in 
Jacobi coordinates [3l] , its numerical solution in SMEC has not yet been given. Firstly, 
if the residual interaction is a contact force, then the interaction of two protons in the 
continuum leads to the ultra-violet divergence. In principle, this divergence could be 
avoided using a finite-range residual interaction. However, such an interaction leads to 
non-local channel-coupling potentials and, hence, equations for channel-projected ujj'^J 
functions become integro-differential coupled-channel equations. 

The three-body cluster model: the Hyperspherical Adiabatic Expansion method In this 
approach, the dynamics of [AJ-particle system is reduced to a three-body dynamics which 
is described by solving Faddeev equations |126lll27j . The three-body wave function is a 
sum of three Faddeev components, each of them expressed in one of the three possible 
sets of Jacobi coordinates {xj,yj}. Each component is then expanded in terms of a 
complete set of hyperangular functions. Essential quantities in this approach are the 
effective radial potentials in the hyperradius variable: 



with rfi^ = (rj — r^y, where Tj is the coordinate of particle i, M = I]?=i ^« and m is 
an arbitrary normalization mass. The set of indices k) is a permutation of (1, 2, 3). 
Xj, are proportional to the distance between two particles and the distance between 
their center of mass and the third particle, respectively. 

The behaviour of effective radial potentials (the hyperspherical adiabatic potentials) 
as a function of p for different angular momentum quantum numbers determines the 
structure of a three-body system at both small and large distances. For short-range 
interactions, the lowest hyperspherical adiabatic potential usually dominates in the 
expansion of the wave function. If both short and long-range interactions are present, 
the wave-function expansion at large distances requires more components. With the 
dominating adiabatic potential, the width of a given three-body resonance can be 




i<k 



(17) 



Two-proton radioactivity 



34 



estimated as 

fPO 

exp 




. (18) 

where Er is the energy of the resonance, pi and po in Eq. f|T8l) are the inner and outer 
classical turning points defining endpoints of the pass through the barrier. If crossings 
of various adiabatic potentials occur along the way through the barrier, the system is 
supposed to follow a path corresponding to the smoothest adiabatic potential. 

The Hyperspherical Adiabatic Expansion (HAE) method has been extensively 
used to study different possible situations in three-body decay resulting from various 
combinations of short and long-range interactions and different binding-energy situations 
of all two-body subsystems |134[ 11531 1154j . In particular, it has been applied for the 



description of the Borromean nucleus ^^Ne(^^0 + p + p) |134l [38] . In this application, a 
parameterized proton-proton interaction was used |155j and the proton-core interaction 
was described by an ^-dependent potential with central, spin-spin and spin-orbit radial 
potentials. 

The HAE approach allows to distinguish between sequential and direct three-body 
decays. Large separation between all pairs of particles in the final state can be achieved 
either by a uniform separation of all particles or by first moving one particle to infinity 
while the others remain at essentially the same distance from each other and second by 
moving the two close-lying particles apart. The latter scenario defines the sequential 
process where the second step of a two-step process starts after the first particle is at a 
distance larger than the initial size of the three-body system. In ^^Ne, both the ground 
state (l/2~) and the first excited state (3/2~) have a lower energy than the proton- 
unstable ground state of ^^F. Hence, the standard sequential decay process through 
the intermediate resonance(s) in [A — 1] nucleus is not possible. However, nothing 
prevents the decay from proceeding via the energetically favourable path described by 
the lowest adiabatic potential until at some point the energy conservation dictates that 
also this two-body structure must be broken. This mechanism which resembles the 
virtual sequential decay process discussed above (see Sect. 18.2. ip leads to a significant 
reduction of the 2p partial half-life of the 3/2~ state in ^''Ne |134l [38]. 

Principal assumptions of the HAE approach are (i) the adiabatic motion in the 
three-body system, and (ii) the cluster ansatz for the many-body wave function of the 
total system. The latter assumption will be discussed in Sect. 18.2.3! in connection with 
another three-body cluster model, the so-called Systematic £^ Model (S£^M) [1281 1129j . 
The adiabaticity assumption is better satisfied for certain collective states and for a 
spontaneous fission process. However, this hypothesis becomes questionable if several 
adiabatic potentials are crossing as found in charged-particles/fragments emission |134[ 



The three-body cluster model: the Systematic i'^-Model The S£^M of Grigorenko et 
al. [1281 1129] is a fiexible tool to simulate various consequences of three-body dynamics 
in 2p emission decay. This model is a variant of a three-cluster model with outgoing fiux. 
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Like in HAE approach, the S^^M uses hyperspherical harmonics to treat the decay of 
a three-body system. However, unhke the HAE method, the S£^M includes essentially 
only the direct three-body decay. In general, the three-cluster models are useful for 
understanding the structure of certain states in light nuclei, such as ^He, ^Li, ^Be |156j . 
^Li, [1571 [158], and ^^Li [156] . 

The basic assumption in S£^M is that the intrinsic constituent particle degrees 
of freedom can be implicitly treated through the judicious choice of both the active 
few-body degrees of freedom and the appropriate effective two-body interactions. 
The cluster ansatz is physically more relevant for states close to the cluster-decay 
threshold [1411 l32l 11421 11431 1144] . Hence, the three-cluster models are perhaps better 
suited for halo systems. On the other hand, couplings to the decay channels involving 
intrinsic states of a daughter nucleus with excitation energies up to ~ 12 MeV [159] 
are strongly present for near-threshold states. Therefore, even in halo states the wave 
function is a superposition of several different cluster configurations. 

In heavier nuclei, the halo wave functions are too complex to be reduced into a 
simple three-cluster parameterization. Nevertheless, the S^^M provides in those nuclei 
an order of magnitude estimate of the relation between the three-body decay energy, 
the effective interactions in the proton-proton and proton-core subsystems, and the 
probabilities of the i"^ configuration of the 2p emitter. Such estimates are very helpful in 
phenomenological applications. Moreover, by solving decay kinematics in hyperspherical 
variables, the S£^M allows to study angular correlations and energy/momentum 
correlations between decay products which presently cannot be investigated by any 
other approach. 

In S£^M, one begins by calculating the internal wave function of the system \E'int at 
a discrete value of the energy E^^t. This wave function is a solution of the Schrodinger 
equation: 

(i^-^i„t)^int = 

with zero boundary condition at some finite value of the hyperradius pint- The time- 
independent part of the wave function for decay particles \E'*-+-' with purely outgoing 
asymptotics and an energy Ei^t is found by solving the inhomogeneous equation 

(i^-Ei,o^^+^ = -^r/2^int 

with the source term given by an internal wave function. This way of modelling 
is based on the assumption that the decaying state is a narrow resonance with a 
negligible coupling to one- and two-particle continua of decay channels. This assertion 
is questionable for couplings to one-particle decay channels, either closed or opened, 
which induce the dynamical effects in the 2p decay. 

The 2p separation energy in S^^M is adjusted using a weak short-range three-body 
force. It is assumed that this force does not modify the barrier penetration process. An 
essential element of the Si'^M is the proton-proton final-state interaction which is tuned 
to the proton-proton phase shift. In the proton-core channel, adjustable potentials of 
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the Woods-Saxon form are employed. The internal structure of the parent nucleus [A] 
is tuned by changing parameters of the proton-core interaction potential for different 
i values. In this way, one generates wave functions with different populations of 
configurations. Clearly, this treatment excludes genuine many-body effects, such as 
the configuration interaction and configuration mixing, or the Pauli blocking which is 
approximated in S^^M by repulsive cores for occupied orbitals. 

Both S£^M and HAE approaches are precisely tuned to describe the three-body 
structures and decays. This gives them a certain conceptual advantage over SM- 
like approaches. On the other hand, an internal structure of the decaying state 
is oversimplified, the treatment of Pauli's exclusion principle is incomplete, and the 
information about the complexity of many-body resonance wave function is lost in 
adjusted interaction potentials. 

Certainly, some of deficiencies of three-body cluster models can be hidden by a 
judicious choice of phenomenological potentials. However, the calculated results can 
hardly be more reliable than the basic assumptions. Due to an unrealistic structure of 
the decaying 2p pair in S^^M, the relative weight of various ^ -configurations obtained 
by fitting the experimental 2p half-life cannot be used to extract meaningful information 
about proton-proton correlations inside a nucleus. It is then astonishing that the S£^M 
can give an order of magnitude estimate of the experimental data. 

9. Comparison between experiment and model predictions 

Results of configuration-interaction approaches (SMEC and extended R-matrix 
approach) can be compared to evaluate the role of nuclear dynamics and continuum 
couplings in the 2p decay. Unfortunately, no meaningful comparison is possible between 
SM-like approaches and three-body cluster models (HAE approach or S£^M). 

In this section, we shall discuss results of SMEC and extended R-matrix approach 
for ^^Fe, ^^Ni and ^^Zn assuming the diproton scenario of 2p radioactivity. We shall also 
evaluate in SMEC an importance of the virtual sequential process in these decays. The 
experimental evidence for the 2p radioactivity in ^^Fe, ^^Ni and ^^Zn has been discussed 
in Sects. 6.2, 6.4 and 6.3, respectively (see also Table 1). 

9.1. The decay o/^^Fe 

In SMEC studies of ^^Fe, internal many-body states in Q are obtained by solving the 
standard SM with the lOKIN Hamiltonian |160] . Configurations with excitations 
from the sd shell to the fp shell are excluded. The residual couplings between Q 
subspace and the embedding continuum is given by the Wigner-Bartlett interaction: 
yires) ^ y^^^ ^ /3P'"]5(r2 - ri), with P = 0.27 {a + p = 1) and Vq = -900 MeV-fm^. 
P°" is the spin-exchange operator. 

The calculated diproton half-life is IS.Sl^gms [31]. The error bar is associated 
with an experimental uncertainty on the decay energy Q2p. External mixing of SM 
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wave functions due to the Q-T couplings gives a negligible correction. More important 
are the Q-V couplings. The external mixing generated by these couplings reduces the 
diproton half-hves by about 10% for Qip in the interval — O.lMeV < Qip < O.lMeV. This 
reduction depends weakly both on the magnitude of Qip and its sign. The experimental 
partial half-life for 2p radioactivity is by factor ~ 4 shorter (cf Table 1). 

The extended R-matrix model yields the diproton half-life of 46lxgms [29]. The 
nuclear structure enters in this model through ^He spectroscopic factors which are 
calculated in fp shell using the GPFXl Hamiltonian [1611 1162] . The calculation using 
FPD6 interaction |163] yields a similar value of the half-life. This half-life is by a factor 
~ 4 slower than the SMEC prediction. 

The sequential 2p decay is always a part of the 2p decay, independently of the 
sign and magnituude of Qip. One should remind that the formal separation between 
the diproton decay and the virtual sequential emission cannot be justified if these two 
modes yield comparable partial decay widths, i.e. T^2p^^ — '^2p''^- 

Let us consider the sequential emission of ^^Fe through the 2~ and 3~ continua 
associated with the ground state {J'^ = 2^) and the first excited state (J'^ = 3") of 
^''Mn, respectively. For Qip < Q2p, both J'^ = 2f and J'^ = 3i states are resonances 
which decay by Ip emission. Details of calculations can be found in Ref. [3^ . 

10 F ' ^ ' ' ' ^ ' ' ' ' —\ 




8 [MeV] 

Figure 17. The energy distribution of the first emitted proton in a sequential decay of 
= 3/2+ ground state of ^^Fe. SMEC calculations, including effect of an external Q- 
V mixing, have been performed for different energies of J"" = 2" ground state of '^'^Mn: 
Qip = -0.1 MeV (solid line), 0.05 MeV (dashed line) and 0.23 MeV (dashed-dotted 
line). The vertical hne denotes Qip — (32p/2. From Ref. [34] . 
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The sequential emission half-life depends strongly both on the position of Ip 
emission threshold (Qip) and on the total decay energy {Q2p)- Assuming that ^^Fe is 
stable with respect to Ip emission {Qip = — O.lMeV), the virtual sequential 2p half-life 
is 165158 about 10 times slower than the diproton decay. This estimate includes 

an effect of external mixing of J'^ = 3/2"'" SM wave functions, which reduces the virtual 
sequential decay half-life by about 30%. For Qip = -|-0.05MeV, even though the standard 
sequential process through an intermediate resonance is possible, its contribution is 
totally screened by the Coulomb barrier. The dominant decay mechanism remains the 
virtual sequential process through the 'ghost' of ground state of ^^Mn, far away from the 
resonance region Qip ± r2-/2 (see Fig. [TTj). The half-life T1/2 = 110^*^29 is about a 
factor of 8 slower than the diproton decay. Couplings to the decay channels (3~, Ij)^^'^^ 
associated with an excited state of '^^Mn are relatively important and further reduce 
the half-hfe by about 30% for Qip = 0.05 MeV and about 38% for Qip = -0.1 MeV. 
One cannot exclude that also decay channels connected with higher lying states of ^^Mn 
contribute significantly to the reduction of the 2p partial half-life |159] . 




- 100 100 200 300 400 500 600 700 800 

Gip [keV] 

Figure 18. The half-lifcs for sequential 2p decay (full points) and Ip decay (squares) 
of ■*^Fe. SMEC calculations include effects of an external Q-V mixing. The dashed- 
dotted line shows the half- life for a diproton decay. From Ref. [34] . 

An energy of the first emitted proton is limited by Q2p. The most probable energy 
of first proton in the virtual sequential decay [3l] {ei^p"^^^ < 0.55 MeV) slightly decreases 
with increasing Qip (see Fig. [TTj). A small asymmetry seen in this distribution is a result 
of different barriers for proton emissions in ^^Fe and in ^^Mn. The full width at half 
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maximum of a ghost is T^ahost) ^ g ^97 MeV at Qip = -0.1 MeV. T^ahost) 

increases with 

Qip in the range of Qip values where the virtual 2p decay dominates. Experimental 
observation of symmetric energy distribution of individual protons [110] indicates that 
either Qip in ^^Fe is less than ~ —0.25 MeV, or the diproton decay mechanism 
dominates. The latter scenario is favoured for typical effective SM interactions and 
Q-values. One should stress that these conclusions have been obtained assuming a 
coherent contribution of diproton and virtual sequential decay mechanisms. 

For Qip < 0.2 MeV, the sequential 2p decay is predominantly related to the strength 
of a ghost in the energy interval ejf - r(^''^''^*)/2 <e < ejf + F^^^"^*) /2 and not to 
the strength of Ip resonance at Qip. In this regime, the sequential 2p decay width does 
not reduce to the product of the width for the first step (Fip) and the branching ratio 
for the second step, as one would obtain semi-classically. Above Qip ~ 0.23 MeV, the 
sequential decay is dominated by the standard two-step process via the ground state of 
^^Mn. 

The Qip- dependence of the sequential 2p decay half-life is shown as full points in 
Fig. [m Squares in Fig. [18] represent the Ip decay half-life and the diproton half- 
life is depicted by the dashed-dotted line. The sequential 2p emission and the diproton 
emission, considered as independent, non- interfering decay paths, yield comparable half- 
lives for 0.15 MeV < Qip < 0.2 MeV. In this interval, the three-body 2p decay cannot 
be reduced to any sequence of two-body decays. 

An interval of Qip values for which the half-lifes of virtual sequential and diproton 
decays are expected to become similar hes close to the range of predicted Qip- 
values [ini [TTl [16]. Hence, both processes could interfere strongly and, in spite of 
symmetric energy distribution of individual protons [llOj . the measurement of Qip 
in ^^Fe is mandatory before drawing definitive conclusions about the 2p emission 
mechanism in this nucleus. 

The sequential 2p decay has two distinct dependences of the decay half-life on Qip 
(see Fig. [18]). For small positive values of Qip (0 < Qip < 0.2 MeV) as well as for 
closed Ip decay channels, the sequential 2p decay half-life changes linearly with Qip. 
This is the regime of virtual sequential 2p decay, where T2p'^^ ^ Fip and the 2p decay 
goes via the 'ghost' of the ground state. For Qip > 0.23 MeV, one enters in the regime 
Fip ^ ^2'p'^\ where the role of 2j~ resonance in "^^Mn is dominant. In this regime of 
standard sequential decay, the 2p partial half-life: 

Ti/2{Qip) ~ exp -(Qip - Qip )^ 
has a minimum at Qip ~ Q2p/'^- 
9.2. The decay o/^^Ni 

In Ref. [31], several SM Hamiltonians including lOKIN [I60], KB3 [I6l| or GXPFl ^THTl 
I162j were used to evaluate the sensitivity of SMEC results to the choice of an effective 
interaction. For the lOKIN Hamiltonian, the diproton half-life is 6.2l2i5nis. The error 
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bars are related with an experimental uncertainty of the decay energy Q2p. External 
mixing of SM wave functions generated by the Q-T and Q-V couplings are unimportant, 
reducing the diproton half-life by ~ 2% in the interval — O.lMeV < Qip < O.lMeV. 
Somewhat stronger is the dependence of 2p partial half-lifes on the SM interaction. For 
GXFPl and KB3 interactions, one finds 7.4l2!9iiis and 6.9t2;7iiis, respectively. For all 
three SM interactions, the SMEC results are consistent with the experimental partial 
half-life. The extended R-matrix model with GXFPl interaction yields a diproton half- 
life 16t4''ms [22] which is about a factor of 2 slower than the SMEC result. 

In the description of virtual sequential emission, the 2p emission through the 
3/2~ and 7/2~ continua associated with the ground state {J'^ = 3/2~) and the first 
excited state (J'" = 7/2~) of "^^Co were considered. SMEC calculations were performed 
with the lOKIN Hamiltonian. Principal uncertainties of the SMEC calculations are 
(i) an unknown position of Ip emission threshold in ^^Ni, and (ii) an insufficient 
precision of the total decay energy. If Qip = —O.lMeV, then the virtual sequential 
2p half-life is 25A~^l'^Q'^ms, i.e. about 4 times slower than the diproton decay. This 
estimate includes an external mixing of 3/2"'" SM wave functions generated by the Q- 
V couplings. These couplings reduce the virtual sequential decay half-life by about 
15%. For Qip = +0.05MeV, the standard sequential decay through an intermediate 
resonance in ^''Co is totally screened by the Coulomb barrier. Hence, the dominant 
decay mechanism remains the virtual sequential process through the 'ghost' of ^^Co 
ground state. The half-life in this case (T1/2 = 16.5l5;7ms) is only ~ 2.5 slower than 
the diproton decay. Couplings to the decay channels (7/2~,/j)°^ associated with the 
excited state of '^''Co are relatively unimportant, reducing the half-life by ~ 6% - 9%. 

9.3. The decay o/^^Zn 

The diproton half-hfe of ^^Zn in SMEC with GXFPl Hamiltonian (TEH [I62] is 
13.81? ;^ms. The error bar is due to the experimental uncertainty on the decay energy 
Q2p iMj- With KB3 Hamiltonian, the diproton half-life is somewhat longer and equals 
17lg°3^ms |164j . The experimental partial half-life for 2p radioactivity is about a factor 
of 4 faster. As in ^^Fe and "^^Ni, a significant contribution from virtual sequential 2p 
decay cannot be excluded but, again, the quantitative estimate of this contribution 
cannot be made without knowing Qip in ^^Zn. 

The extended R-matrix model with GXFPl interaction yields the diproton half-life 
1014ms j2l] which is about 30% shorter than the SMEC prediction. 

9.4- Discussion 

The coupling of a nuclear system to the environment of its decay channels may change its 
properties in a non-perturbative way |32]. The magnitude of the system - environment 
coupling depends explicitly on the location of various emission thresholds and on the 
structure of 5'-matrix polers. The description of resonances and their decays are mainly 
sensitive to the two ingredients: the effective Hamiltonian TCqq (see Sect. 18. ip and 
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the unitarity of the S'-matrix which causes a non-trivial energy dependence of the 
couphng matrix elements between resonance states and continuum. TYqq includes basic 
structural information contained in Hqq (the closed quantum system Hamiltonian) and 
the coupling matrix elements between discrete and continuous states which depend on 
the structure of the decay channels. 

SMEC studies of 2p decays in "^^Fe, "^^Ni and ^^Zn have shown that the virtual 
sequential 2p emission cannot be separated easily from the diproton decay, even though 
they describe two different limits. The first one corresponds to two uncorrelated proton 
emissions separated in time, whereas the latter one represents an instantaneous emission 
of strongly correlated two protons. The most probable energies of these protons are 
nearly the same in both scenarios. 

If half-lives for these two pathways are strongly different, then neglecting their 
mutual interference and the omission of spatio-temporal correlations in the virtual 
sequential decay can be probably justified. Otherwise, the 2p emission described by 
an effective Hamiltonian (JTTi) has to be considered as a genuine three-body decay. 
Unfortunately, the Ip emission thresholds in ^^Fe, ^^Ni and ^"^Zn are unknown so that one 
cannot rule out the possibility of a strong interference between (virtual) sequential and 
diproton decay modes. SMEC calculations in "^^Fe with Qip values which are consistent 
with available theoretical estimates [151 [13 IIS] show that the virtual 2p sequential 
half-life is by a factor of about 10 longer than the diproton half-life. This would favour 
an interpretation of the observed 2p emission in this nucleus as predominantly the 
diproton decay. Such a conclusion agrees also with the observation of a symmetric 
energy distribution of individual protons in the decay of '^^Fe |110j . 

For the sake of argument, let us suppose that the virtual sequential emission 
mechanism yields much longer half-lives in all three nuclei: ^^Fe, ^^Ni and ^^Zn. It 
is then interesting to notice that the diproton decay mechanism in SMEC reproduces 
quite well the experimental partial 2p decay half-life in a doubly magic nucleus (^^Ni) 
and fails by a factor of about 4 in both open-shell nuclei (^^Fe and ^^Zn). If this 
tendency is not coincidental, or follows from a poor knowledge of Q-values and an 
insufficient experimental statistics, then one might be tempted to associate it with 
pairing correlations. As suggested in Ref. [21], the change of pairing correlations 
induced by the configuration mixing beyond the IpOf shell may well be responsible 
for this enhancement. A similar effect has been seen a long time ago in (p,t) transfer 
reactions [165] . The absence of this component of pairing correlations is more important 
in open-shell nuclei than in a closed-shell nucleus "^^Ni. Hence, one may expect that the 
fp shell SMEC calculations are better suited to describe the observed 2p partial half-life 
in ^^Ni than in "^^Fe and ^^Zn. 

One should mention that different effective interactions in the fp and sdfp shells 
give similar results to within ~ 10 — 20%. External mixing, generated by the Q-V 
coupling, modifies the diproton decay width by less than 10% in ^^Fe and ~ 2% in ^®Ni. 
Again, the SMEC predictions for the closed-shell nucleus ^^Ni seem to be more reliable 
since the direct Qip-dependence of the diproton half-life is weaker in this case. 
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The R-matrix model [2S] calculations show a gradual decrease of the 2p partial 
half-life going from ^^Fe to ^^Zn across the shell closure at Z=28, N=20. Consequently, 
the diproton half-life in '^^Fe is by a factor of ~ 5 longer than in ^^Zn. On the contrary, 
both in the SMEC results and in the data, the 2p partial half-lives of ^^Fe and ^'^Zn are 
very close. 
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Figure 19. The partial lialf-livcs for the diproton decay from the ground state of 
2gFe, 2i-Ni, and |oZn, as calculated in the generalized R-matrix theory [29] (full 
triangles) and in the SMEC [33] (open circles). The experimental values are shown 
with black points. For reasons of clarity, the uncertainty of theoretical results due to 
the uncertainty of the Q2p values and experimental error bars are not shown. 

Fig. [T9] compares the experimental partial half-lives for the diproton decay of 26^^) 
2gNi, and IqZu with the calculated ones in the SMEC and in the extended R-matrix 
model. One may notice two characteristic features from this comparison. Firstly, SMEC 
and R-matrix models predict a different Z-dependence of calculated 2p partial half-lives 
around the doubly-magic nucleus ^^Ni. Secondly, SMEC yields an opposite tendency of 
^^Fe and ^^Zn half-lives with respect to the half-life of ^^Ni than seen experimentally. 
Concerning the first point, one may tentatively assume that this feature is related to 
the continuum coupling induced dynamical effects in the 2p decay which are missing in 
the R-matrix description. Certainly, one does not expect that this different systematic 
tendency of calculated half-lives in the SMEC approach and in the R-matrix model 
could be explained by the absence of configuration mixings beyond the IpOf shell, since 
both in SMEC and in R-matrix calculations these configurations are not included. 
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Concerning the second point, it may be related to the small number of channels in 
SMEC calculations of the Q-P-coupling induced an external configuration mixing. In 
Ref . [3l] , the external mixing is a result of coupling between the ground state of nucleus 
[A] and the two channels corresponding to the ground state and the first excited state of 
a nucleus [A — 1]. This restriction on the number of channels is justified in closed-shell 
nuclei as in this case the ground-state to ground-state coupling of the [A] and [A — 1] 
nuclei exhausts up to about 90% of the total coupling strength. For open-shell nuclei, 
the strength of continuum coupling is largely spread over many channels, both open 
and closed. Indeed, the SMEC studies of continuum coupling effects in oxygen and 
fluorine isotopes |159j have shown that a contribution of the ground-state to ground- 
state coupling to the total continuum coupling strength may vary between a few % 
and ~ 30% in open-shell nuclei. The size of a model space involved in the calculation 
of the 2p decay in 26^^, 28^1) ^ind lo^n makes it impossible to verify this observation 
numerically, but if true it may lead to a strong enhancement of the contribution from 
the virtual sequential 2p decay path and, henceforth, to a decrease of the 2p partial 
half-life. One should remind, that the interference between the virtual sequential and 
diproton emission paths cannot be probably neglected if both processes become equally 
important. 

As said before, SMEC predictions are very sensitive to Q values which for the nuclei 
studied are either not known experimentally with a sufficient precision {Q2p values), or 
are unknown {Qip values). The virtual sequential decay half-lives depend explicitly and 
strongly on the continuum coupling dynamics between the Q, V and P, T subspaces 
which is governed by Q values. A direct dependence of the diproton half-lives on Q-V 
virtual couplings is rather weak, i.e. dressing of the diproton decay mode by virtual 
Ip excitations does not change significantly its properties in ^^Fe, ^^Ni and ^^Zn 0. On 
the other hand, this idealized decay mode represents a sensible approximation of the 
2p decay process described by (IHl) only if the virtual sequential decay path is strongly 
suppressed. This in turn depends directly on the dynamics of the Q, V and V, T 
virtual couplings. Hence, further progress in describing spontaneous 2p radioactivity and 
understanding the role of both pairing correlations and virtual excitations to continuum 
states in this process depends on a precise experimental determination of the Qip and 
Q2p values. 

10. Conclusions and outlook 

The study of 2p emission opens a new window to investigate nucleon-nucleon correlations 
and the structure of atomic nuclei. Although these studies, at least partially, could 
have been performed already with /3-delayed 2p emitters, the field received a strong 

I This is in contrast to the diproton decay from the second excited 1~ state of ^*Ne In this 

case, the diproton mode is dressed by virtual excitations to the opened Ip emission channels leading to 
the particle stable 5/2^, 1/2^ final states in ^^F. In general, couplings to open Ip emission channels 
strongly modify properties of the diproton mode. 
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boost with the observation of ground-state 2p radioactivity. This discovery triggered 
new experimental developments, such as time-projection chambers, and new theoretical 
activity related mainly to the unified description of structure and reactions with two- 
particle continua and the development of various methods to deal with three-body 
asymptotic states of charged particles/fragments. 

Both experiment and theory have reached a certain level of sophistication, but 
significant improvements are still needed to develop 2p radioactivity into a powerful 
tool for nuclear structure studies. On the theoretical side, the lack of three-body 
asymptotic channels hampers the calculation of various correlation functions in the 
SMEC framework. Description of direct 2p emission with three-body asymptotics in 
SMEC requires inclusion of finite-range interactions. At the moment, the angular 
correlations between decay products can be studied solely in the three-body cluster 
models. The simplified description of nuclear structure in cluster models does not 
allow to draw deeper conclusions about 2p in-medium correlations from these angular 
correlations. 

As one deals with three particles in the final state, the 2p emission can occur 
in different ways, i.e. the asymptotic channel cannot be defined unambiguously. 
Consequently and in contrast to one-proton or a radioactivity, the 2p radioactivity does 
not have a unique experimental signature and its interpretation is always depending on 
the theoretical interpretation. A complete theoretical description of 2p emission has 
to take into account all different decay paths and include interferences between them. 
It is not obvious whether different 2p decay paths, such as the diproton path and the 
virtual sequential 2p decay path which both yield similar proton energy distributions, 
can be distinguished by, e.g., the proton-proton and/or proton- fragment correlations. 
Further studies along this fine are necessary for a best possible characterization of 2p 
radioactivity. 

On the experimental side, more cases of 2p radioactivity have to be studied and the 

quality of the data has to be improved, as the possibilities to compare the experimental 
data existing today with modern theories of 2p radioactivity are hampered by large 
experimental error bars. To improve the experimental data, increased statistics is only 
one aspect. In the case of ^"''Fc, the major contribution to the uncertainty of the decay 
energy comes from the energy calibration. Therefore, new and improved calibration 
tools have to be developed as well. A precise determination of Qip and Q2p would 
help to check different features of SMEC calculations, such as the choice of effective 
interactions and valence space, or the absence of explicit three-body asymptotics. 

All experiments performed up to now seem to indicate that the sequential decay 
prevails, as long as intermediate states for Ip emission are energetically accessible. The 
only exceptions seem to be the decays of excited states in ^''Ne and in ^^Ag"^, where 
despite many intermediate states a strong proton-proton correlation was observed. A 
halo-like structure and/or dynamical effects associated with the tunnelling through the 
anisotropic Coulomb barrier could be at the origin of the observed correlation. Higher 
statistics data and advanced theoretical studies are needed to confirm or reject this 
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interpretation. 

Another application of the 2p studies are direct resonant 2p radiative capture 
processes |166l 1167] . It was demonstrated that the astrophysical 2p capture rates can 
be enhanced by as much as several orders of magnitude in certain ranges of excitation 
energies. Similar tendency was found for the size of a nonresonant El contribution to 
the 2p capture rates |167j . These predictions may have important consequences for the 
rp-process path changing earlier predictions concerning rp-process waiting points [168j . 

The advent of new radioactive beam machines like the RIBF facility in RIKEN, 
Japan or the FAIR facility in GSI, Germany should give access to new ground-state 
2p emitters such as ^^Ge, ^^Se, or ^^Kr and allow one to achieve much higher rates 
for the known ones. These possibilities will most likely significantly increase the body 
of experimental data and may establish 2p radioactivity as a powerful tool of nuclear 
structure studies. 

These new facilities may allow also to search for a similar phenomenon which is 
one- and two-neutron radioactivity. From an experimental point of view, it seems to 
be rather difficult to find good candidates for these new decay modes. Due to the 
absense of the Coulomb barrier for neutrons, only the centrifugal barrier will prevent 
these neutrons from escaping "immediately" from the nucleus. To create a sufficiently 
high barrier, high-j orbitals are needed which form the nuclear ground state only in 
medium-mass nuclei. However, these medium-mass neutron drip-line nuclei are rather 
far away from the valley of stability and extremely difficult to produce. Possibly, high-j 
isomers could be candidates for this kind of radioactivity. 

From a theoretical point of view, two-neutron emitters would be ideal cases to 
study. All theoretical formalisms can be applied to the two-neutron case without major 
changes. Actually, the two-neutron three-body problem can be solved exactly and not 
approximately as in the three-body case with charged particles. Moreover, angular 
correlations will be directly related to the internal two-neutron correlations without the 
uninteresting correlations due to a possible anisotropy of the Coulomb barrier. 
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